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CHAPTER  1: 
Introduction 


The  process  of  finding  the  location  of  one  or  more  unknown  subjects  is  known  as  Search  and 
Detection  Theory  (SDT)  and  is  of  great  interest  to  U.S.  Navy  (USN).  While  there  are  many 
applications,  both  military  and  non-military,  SDT  is  of  particular  import  to  Anti-Submarine 
Warfare  (ASW).  This  is  due  to  the  desire  to  locate  adversarial  submarines  because  the  nature 
of  submarines  is  to  use  their  covert  operating  patterns  to  conduct  Intelligence,  Surveillance 
and  Reconnaissance  (ISR),  Anti-Surface  Warfare  (ASUW),  Strike  Warfare  (SW)  and  Special 
Operation  Forces  (SOF)  missions. 

With  the  advent  of  high  end  diesel-electric  submarines,  increased  tensions  and  activity 
from  capable  nuclear  submarines  (Russia  and  China),  and  the  proliferation  of  submarine 
technologies  to  previously  non-capable  navies,  the  need  for  better  approaches  to  SDT  is 
evident.  The  proliferation  of  these  submarine  technologies,  and  the  rise  of  autonomous 
vehicles  such  as  unmanned  underwater  vehicles  (UUV)  and  unmanned  surface  vehicles 
(USV),  both  as  preprogrammed  searchers  and  potential  targets,  dictate  that  the  parameters 
of  SDT  algorithms  must  be  investigated. 


1.1  History 

The  Operations  Evaluation  Group  (now  known  as  the  Center  for  Naval  Analyses  (CNA))  [1] 
was  created  during  World  War  II  (WWII)  to  enlist  the  help  of  mathematicians  and  scientists 
to  assist  in  developing  new  ASW  tactics  and  strategies.  This  small  group  of  MIT  scientists 
worked  with  ASW  Operations  Research  Group  (ASWORG)  directly  in  the  field  for  optimum 
observations  and  real-time  solutions  [2].  These  scientists  used  the  burgeoning  search  theory 
to  improve  U-boat  search  and  neutralization,  especially  with  favorable  patterns  and  operating 
characteristics  for  patrol  aircraft. 

After  the  war,  search  theory  was  extended  to  other  naval  operations  such  as  searching 
for  the  USS  Scorpion  (SSN-589),  lost  ordnance,  and  tracking  Russian  ballistic  missile 
submarines  [3].  The  United  States  Coast  Guard  (USCG)  developed  a  computer-aided 
search  program  (CASP)  [4]  for  missing  boats  and  mariners  lost  at  sea. 
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Dr.  Bernard  Koopman  of  the  WWII  ASWORG  developed  the  historical  basis  for  search 
theory  by  dictating  the  following  basic  principles. 

1.  Quantify  the  probability  of  the  sensor  capability  to  detect  the  target. 

2.  Classify  the  movement  of  the  searchers  and  target. 

3.  Analyze  the  constraints  on  the  searcher’s  forces  and  resources 

4.  Develop  the  highest  probability  of  detection  with  the  greatest  economy  of  forces  [5]. 

1.2  Applications  and  Previous  Research 

The  applications  for  SDT  are  numerous,  so  that  while  the  context  of  this  research  is  limited 
to  an  ASW  submersible  search  problem,  the  basics  can  be  extended  to  various  problem  sets 
and  models.  Research  conducted  at  Naval  Postgraduate  School  (NPS)  has  focused  on  naval, 
non-naval,  and  other  Joint  Warfare  (JW)  applications.  Research  has  also  been  conducted  at 
non-military  institutions  for  civilian  applications. 

1.2.1  Some  Previous  Research 

1.  Zachary  Lukens’s  2016  thesis  [6]  investigated  SDT  applied  to  a  single  target  in  the 
South  China  Sea  and  two  searchers  using  cookie  cutter  sensors  utilizing  inline  spiral, 
inline  ladder,  and  multipath  ladder  searches.  His  research  found  that  increasing  either 
the  detection  radii  or  speeds  decreased  the  time  to  detection  and  that  spiral  search 
was  faster  than  ladder  searches  at  large  radii,  but  both  ladder  searches  were  faster  than 
spiral  search  for  small  radii. 

2.  Volkan  Sozen’s  2014  thesis  [7]  used  SDT  applied  to  Unmanned  Aerial  Vehicles  (UAV) 
patrolling  a  border.  He  modeled  the  search  as  one  UAV  searching  the  area,  two  UAVs 
searching  the  same  area,  and  two  UAVs  independently  searching  two  equally  split 
sub-areas.  His  research  found  that  detection  increased  with  the  second  searcher 
and  splitting  the  searchers  into  independent  sub-areas,  the  border  was  covered  more 
efficiently. 

3.  Dr.  Michael  Atkinson  [8]  (together  with  Moshe  Kress  and  Rutger- Jan  Lange)  re¬ 
searched  searching  for  a  target  discrete  location  based  on  receiving  of  intelligence 
of  varying  reliability  and  the  time  cost  incurred  based  on  searching  a  certain  area. 
His  research  became  the  problem  of  looking  at  the  probability  of  detection  based  on 
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amount  of  intelligence  gathered  weighed  against  the  target  escaping  due  to  a  delayed 
search.  This  has  direct  applications  to  military  and  non-military  problems  based 
on  the  information  quality.  The  study  concluded  that  high  intelligence  credibility 
had  better  results  in  waiting  to  attack,  while  low  credibility  had  better  results  for 
immediate  attack. 

4.  H.  Zhou  and  H.  Wang  in  [9]  examined  optimum  paths  for  a  solitary  searcher  looking 
for  a  solitary  target  using  spiral,  ladder,  loop  and  square  spiral  paths. They  showed 
that  a  stationary  searcher  is  optimally  placed  in  the  center  of  the  search  area.  It  also 
showed  that,  with  a  moving  searcher,  the  ladder  pattern  gives  the  most  likely  detection. 
In  [10],  they  examined  optimum  searcher  placement  seeking  a  target  travelling  back 
and  forth  between  a  detectable  operating  area  and  an  undetectable  hiding  area.  They 
showed  closed  form  solutions  for  determining  mean  time  to  detection  based  on  the 
placement  of  the  searcher  in  a  given  path  from  hiding  to  operating  area  or  within 
the  operating  area  itself.  Their  results  showed  how  to  distribute  sensors  ("searchers") 
among  the  routes  and  the  operating  area.  They  extended  their  study  to  include  sensors 
in  detecting  a  target  moving  between  a  hiding  area  and  an  operating  area  [1 1],  Based 
on  the  average  time  to  detection,  they  were  able  to  evaluate  the  performance  of  placing 
the  sensor(s)  to  monitor  various  travel  pathway(s)  or  to  scan  the  operating  area. 

1.2.2  Applications 

The  US N  has  utilized  SDT  for 

•  ASW 

•  ASUW 

•  Anti-Air  Warfare  (AAW) 

•  SW 

•  Mine  Warfare  (MIW) 

•  other  operations. 

Other  branches  and  JW  applications  include 

•  ISR 

•  Counter-drug  operations 

•  Counter-terrorism  operations 


3 


•  Ballistic  Missile  Defense 

•  SW  and  Artillery  Fires  operations 

•  Search  and  Rescue. 

Multiple  civilian  applications  are  also  prevalent  in  research  and  development,  to  include: 

•  Search  and  Rescue 

•  Quality  Control 

•  Resource  exploration  and  search  (oil,  minerals,  fishing,  etc.) 

•  Microeconomics  (resource  to  market  matching). 
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CHAPTER  2: 
Discrete  Search 


2.1  Theory 

The  following  sections  outline  the  basic  concepts  and  theory  behind  the  discrete  search 
theory. 

2.1.1  Terms 

A  -  Total  search  area  (no  searches  conducted  outside  of  A) 
n  -  Total  number  of  equal  area  subdivisions  (boxes)  of  A  :  (i  £  {1../?}) 

Pi  -  /J[tgt  in  / th  box]  probability  of  target  present  in  /th  box 

qi  -  P[miss  detection  |tgt  in  /th  box]  probability  of  failing  to  detect  target  in  /th  box 
given  it  is  present  in  /th  box 

ki  -  number  of  glimpses  (discrete  time  interval)  with  search  was  conducted  in  /th  box 

n 

m  -  total  number  of  glimpses:  m  =  Y  ki 

i—  1 
n 

By  the  concept  of  total  probability,  Y  Pi  =  1. 

1  =  1 


2.1.2  General  Concepts 

The  basic  assumptions  for  the  discrete  search  model  are  [12]: 

•  Single  stationary  target 

•  Single  searcher  with  imperfect  sensor  (<?;) 

•  Total  search  time  limited  (m  glimpses) 

•  No  time  lost  transitioning  from  one  box  to  another 

•  Each  search  conducted  independent  of  other  glimpses 

•  No  false  positives 

•  Non-adaptive  search. 

The  probability  of  the  target  not  being  detected  in  the  /th  box  on  a  given  glimpse  is 
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P[miss  detection]  =  P[miss  detection  |  tgt  in  /th  box]P[tgt  in  / th  box]  =  PlM  =  p,qi  and 
the  detection  probability  on  a  given  glimpse  for  the  /th  box  is  ^[detection]  =  Pfdetection  | 
tgt  in  /th  box]P[tgt  in  /th  box]  =  PlD  =  p,(  1  -  q{).  After  kj  glimpses,  the  probabilities  of 
missing  or  detecting  in  the  /th  box  are  Pff  =  piqkl  and  Pl'ki  =  p,( 1  -  qk‘),  respectively.  This 
leads  to  [13],  by  the  total  probability  concept,  the  likelihood  of  total  failure  and  detection  as 

n 

IJM  =  ^  Piqf‘  (2.1) 

1  =  1 

and 

n 

PD  =  J]pi(l-qfi).  (2.2) 

i= 1 


2.1.3  Model  Overview 

The  greedy  algorithm  is  a  myopic  search  model  that  is  designed  to  attempt  to  optimize 
Equation  2.1  in  order  to  streamline  the  detection  based  on  a  limited  number  of  resources 
(time,  searchers,  etc.).  For  each  glimpse,  the  probability  of  detection  is  calculated  given 
previous  glimpses  have  failed  to  detect  the  target.  This  change  in  detection  probability  [13] 

A  =  (1  -  q-i+1)Pi  -  (1  -  qf)Pi  =  (1  -  qdq^Pi  (2.3) 

is  used  to  select  the  maximum  detection  probability  and  the  search  is  conducted  in  that  box, 
increasing  kj  by  1  for  that  box.  This  updates  the  probability  for  that  box  alone.  If  the  target 
is  not  found,  the  search  proceeds  to  the  next  glimpse  using  the  same  formula.  This  will  be 

n 

repeated  until  the  target  is  found  or  £  kj  =  m.  The  cumulative  detection  probability  is  the 

i=i 

sum  of  the  maximum  A  for  at  every  step  or  Equation  2.3. 

2.1.4  Example 

The  example  in  Table  2.1  from  [12],  outlines  the  "greedy"  search  algorithm.  The  target 
is  in  one  of  three  boxes,  with  associated  location  and  conditional  miss  probabilities.  The 
probability  of  detection  for  each  cell  is  calculated  based  on  finding  the  target  on  the  first 
glimpse  with  no  previous  looks  (k  =  (0, 0, 0))  so  that  qkl  =  1 .  The  highest  probability  is 
in  the  first  cell  making  the  qki  <  1  reducing  the  probability  on  the  next  round  of  glimpses. 
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The  cumulative  probability  after  six  glimpses  is  the  sum  of  the  maximum  probabilities  at 
each  round,  Pd  =  -25  +  .125  +  .1  +  .08  +  .064  +  .0625  =  .6815. 

Table  2.1.  Optimal  Sequence  of  Discrete  Looks.  Adapted  from  [12], 


{k^k^ks) 

/>t(l  -  q\)q\l 

Pi(  i  -  qi)q\2 

/>3(1  -  q*)q\3 

(0,0,0) 

.25 

.08 

.1 

(1,0,0) 

.125 

.08 

.1 

(2,0,0) 

.0625 

.08 

.1 

(2,0,1) 

.0625 

.08 

0 

(2,1,1) 

.0625 

.064 

0 

(2,2,1) 

.0625 

.0512 

0 

Pi 

P2 

P3 

.5 

A 

.1 

q\ 

<72 

<73 

.5 

.8 

0 

2.2  Problem  and  Model 

This  basic  model  for  discrete  search  is  based  solely  on  the  conditional  miss  probabilities 
(qi)  and  the  location  probabilities  (pi).  These  parameters  will  be  varied  and  results  shown. 

2.2.1  Parameter 

The  location  probabilities  will  be  varied  according  to  Figure  2. 1  and  the  detection  probabil¬ 
ities  will  be  varied  according  to  Figure  2.2.  The  uniform  location  and  conditional  detection 
probabilities  are  to  set  up  the  base  case  to  compare  other  distribution  combinations  too. 

The  base  location  probability  per  box,  po  =  1  =  as  shown  in  the  uniform  case.  Two 
linear  stratified  location  probability  densities  are  used  to  simulate  the  target  being  placed 
preferentially  along  some  line  of  longitude  (small  scale)  for  the  first  distribution  or  some 
line  of  constant  longitude/latitude  line.  Each  of  these  increase  by  a  factor  of  10  over  the 
search  area  in  the  direction  of  greatest  change.  Three  Gaussian  distributions  are  used  to 
simulate  the  likelihood  of  the  target  having  a  well-localized  location  probability  prior.  They 

(X-Xp)2  ( IJ-IJQ )2 

follow  a  distribution  function  p(x,  y)  =  — \=e  2a~x  2o>  .  The  centered  Gaussian  has 

cryln 

different  criat  and  <Jioncj.  The  other  two  Gaussian  plots  are  centered  on  a  corner  and  side 
with  criat  =  criong. 
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The  conditional  detection  probability  densities  are  uniform  as  a  base  sensor  that  works 
the  same  in  all  locations  qt  =  .3,  giving  a  70%  of  detection  at  each  location.  The  two 
stratified  detection  probability  densities  are  linear  stratified  along  latitude  or  longitude  from 
q\  =  (.063,  .63).  This  simulates  the  sensor  performance  based  on  some  geographic  profile 
(i.e.,  change  in  bottom  type  for  SONAR). 


Figure  2.1.  Location  Probability  Distributions  pi 


From  top  left  clockwise:  Uniform  d i st . ,  Linear  stratified  along  longitude,  Linear 
stratified  along  Long/Lat,  Gaussian  d i st . ,  Gaussian  centered  in  corner,  Gaussian 
centered  on  side 


Figure  2.2.  Conditional  Miss  Probability  Distributions  qi 


along  latitude 
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2.3  Results 

The  following  sections  show  the  results  for  six  different  location  probability  densities  and 
the  three  conditional  search  probability  densities.  Each  combination  was  examined  for 
n  =  50  glimpses  (half  the  number  of  available  boxes). 

Since  there  are  18  possible  combinations  each  with  various  tables  and  figures,  the  following 
will  show  a  small  sample  of  the  results.  All  of  the  figures  and  tables  will  be  presented  in 
Appendices  A  and  B,  respectively. 

2.3.1  Linear  Stratified  Location 

Figures  2.3,  2.4,  and  2.5  show  the  remaining  location  probability  density  for  the  three 
conditional  search  probability  densities  after  50  glimpses.  Also  shown  are  the  boxes  that 
were  selected  to  search  by  the  greedy  algorithm. 

For  the  uniform  conditional  search  probability  in  Figure  2.3,  unsurprisingly,  the  selected 
boxes  start  with  the  highest  /?,■  and  slowly  work  to  the  left  as  pi  lowers. 


Figure  2.3.  Stratified  Location  and  Uniform  Search  Density  Results 


Remainder  Prob. 

Pi  xl  O-3  Actual  Search  N 

|  =C  O 

6  ^ 

4 

J3  ° 

2 

0.5 

3  0  5 

Long. 

Stratified  a 

Ps 

3  0  5 

XT  rufor  ru 

Cfi 

0.3 

0.2 

O  .  1 

1  ^tj  o 

0.01  5 

001  o 

0.005 

3  0  5 

3  0  5 

Remaining  probability  of  detection  and  boxes  search  based  on  stratified  pt  and 
uniform 


For  the  longitudinal  stratified  conditional  search  probability  in  Figure  2.4,  the  selected 
boxes  also  started  with  highest  p;  and  slowly  work  to  the  left  as  p,-  lowers.  This  is  slightly 
unexpected  since  a  better  probability  of  detection  exists  in  areas  of  lower  target  location 
probability.  But  the  factor  of  1  -  is  93.7%  in  the  western  side  and  37%  in  the  east.  This 
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coupled  with  the  pi  being  a  factor  of  10  higher  in  the  east  leads  it  to  be  a  higher  probability 
of  detection  even  though  the  sensor  performs  worse  there. 


Figure  2.4.  Stratified  Location  and  Longitudinal  Stratified  Search  Results 


Remaining  probability  of  detection  and  boxes  search  based  on  stratified  p, 
and  linear  stratified  qt 


Figure  2.5.  Stratified  Location  and  Corner  Stratified  Search  Density  Results 


Remaining  probability  of  detection  and  boxes  search  based  on  stratified  p,  and 
corner  stratified  qi 


For  the  corner  stratified  conditional  search  probability  in  Figure  2.5,  the  selected  boxes  start 
in  the  southeast  due  to  the  juxtaposition  of  high  pi  and  low  q-,.  The  path  moves  mostly 
along  discretized  radial  arcs  centered  on  the  southeast  corner.  Which  box  is  chosen  via  the 
greedy  algorithm  is  dependent  on  when  1  -  r/,  will  push  a  probability  of  detection  higher 
for  a  lower  p,-. 
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For  the  longitudinal  stratified  p,  location  probability  density,  the  following  results  were 
found  based  on  running  the  search  for  up  to  500  glimpses.  A  data  set  using  a  step  size  of 
5  is  shown  in  Figure  2.6  along  with  associated  fit  lines  as  calculated  from  MATLAB’s  fit 
algorithm  based  on: 

f(x)  =  1  -  e~b\ 

Each  of  the  fit  lines  are  color  coded  to  the  data  set  with  black  for  uniform,  red  for  longitudinal 
stratification,  and  green  for  corner  stratification.  As  the  number  of  glimpses  increases,  the 
probability  approaches  1,  with  the  uniform  case  increasing  faster  than  the  other  two. 

Figure  2.6.  Stratified  Location  Cumulative  Probability  and  Fit  Lines 

Pd  vs  Number  of  Glimpses  and  Exponential  Fitted  Curves 

1 

0.9 
0.8 
0.7 
0.6 
CL?  0.5 
0.4 
0.3 
0.2 
0.1 
O 

O  lOO  200  300  400  500 


Cumulative  probabilities  of  detection  based  on  the  number  of  glimpses.  In¬ 
cluded  are  the  fit  best  fit  line 
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Table  2.2.  Parameter  Variation  for  p,  for  Longitudinal  Stratification 


Parameter 

b 

r2 

SSE 

RMSE 

P  strut  ^uniform 

PstrcitiQlong 

P st rat  ^corner 

0.013208 

0.0089597 

0.011187 

0.99908 

0.99920 

0.99808 

0.00065853 

0.00066179 

0.0013515 

0.008554 

0.0085751 

0.012254 

Best  fit  line  parameters  for  a  longitudinally  stratified  location  probability  density 
of  the  form  fix)  -  I  -  e~bx 


2.3.2  Gaussian  Corner  Location 

The  remaining  location  probability  densities  for  the  three  conditional  search  probability 
densities  after  50  glimpses  are  shown  in  Figures  2.7,  2.8,  and  2.9.  Also  shown  are  the  boxes 
that  were  selected  to  search  by  the  greedy  algorithm. 


For  the  uniform  conditional  search  probability  in  Figure  2.7,  unsurprisingly,  the  selected 
boxes  start  with  highest  pi  and  slowly  work  outward  as  p,  depletes.  A  second  look  begins 
occurring  in  those  higher  probability  boxes  prior  to  searching  all  boxes. 


Figure  2.7.  Gaussian  Corner  Location  and  Uniform  Search  Density  Results 


Remaining  probability  of  detection  and  boxes  search  based  on  Gaussian  p-,  and 
uniform  q{ 


For  the  longitudinal  stratified  conditional  search  probability  in  Figure  2.8,  the  selected 
boxes  also  started  with  highest  p,  and  slowly  works  outward  as  p,  lowers.  This  is  slightly 
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unexpected  since  there  is  better  probability  of  detection  in  some  areas  of  lower  target 
location  probability.  As  the  probability  of  detection  decreases  on  subsequent  looks,  the 
lower  qj  increases  the  likelihood  of  finding  the  target  to  the  west  of  the  initial  location 
locus.  From  the  Remainder  Probability  graph  in  Figure  2.8,  the  higher  location  probability 
remains  dispersed  in  the  northwest  corner  albeit  non-uniformly. 


Figure  2.8.  Gaussian  Corner  and  Longitudinal  Stratified  Search  Density  Re¬ 
sults 
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Remaining  probability  of  detection  and  boxes  search  based  on  stratified  pt  and 
linear  stratified  qt 


For  the  corner  stratified  conditional  search  probability  in  Figure  2.9,  the  selected  boxes  start 
in  the  northeast  due  to  the  high  p,  overriding  the  high  qt.  The  Remainder  Probability  graph 
in  Figure  2.9  shows  some  areas  2-3  times  greater  than  other  areas,  especially  around  the 
original  locus.  Some  areas  will  be  searched  multiple  times  before  every  box  is  searched 
once. 
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Figure  2.9.  Gaussian  Corner  and  Corner  Stratified  Search  Density  Results 


Remaining  probability  of  detection  and  boxes  search  based  on  stratified  p,-  and 
corner  stratified  g,; 


For  the  Gaussian  Corner  p,  location  probability  density,  the  following  results  were  found 
based  on  running  the  search  for  up  to  500  glimpses.  A  data  set  using  a  step  size  of  5 
is  plotted  in  Figure  2.10  along  with  associated  fit  lines  as  calculated  from  MATLAB’s  fit 
algorithm  based  on: 

m  =  1  -  e~b\ 

Each  of  the  fit  lines  are  color  coded  to  the  data  set  with  black  for  uniform,  red  for  longitudinal 
stratification,  and  green  for  corner  stratification.  As  the  number  of  glimpses  increase,  the 
probability  approaches  1,  with  the  uniform  case  increasing  faster  than  the  other  two. 
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Figure  2.10.  Gaussian  Corner  Cumulative  Probability  and  Fit  Lines 


Cumulative  probabilities  of  detection  based  on  the  number  of  glimpses.  In¬ 
cluded  are  the  fit  best  fit  lines 


Table  2.3.  Parameter  Variation  for  pi  for  Corner  Gaussian  Location 


Parameter 

b 

r2 

SSE 

RMSE 

PGC  uniform 
PGCitflong 

PGCiQcorner 

0.03117 

0.01869 

0.01441 

0.99683 

0.9952 

0.99523 

0.0013941 

0.0027014 

0.0029724 

0.012446 

0.017325 

0.018173 

Best  fit  line  parameters  for  a  Gaussian  Corner  location  probability  density  of  the 
form  f(x)  —  I  -  e~bx 


2.3.3  Cumulative  Detection  Probability 

The  total  probability  of  detection  for  the  after  50  glimpses  is  shown  in  Table  2.4.  The 
Gaussian  p,’s  have  higher  probabilities,  which  is  not  unexpected  since  they  have  a  much 
higher  concentration  on  a  few  points.  This  allows  for  multiple  searches  over  a  smaller 
number  of  points  instead  of  50  different  boxes.  The  other  result  seems  to  be  that  higher 
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probability  results  when  p,  matches  up  to  1  -  q\  as  it  does  when  p strati fied  has  a  higher 
probability  of  detection  when  paired  with  qCOrnerstrat  instead  of  qiongstrat. 


Table  2.4.  Pd  after  50  glimpses 


\  ?i 

Pi  "\ 

Uniform  Long.  Stratified 

Corner  Stratified 

Uniform 

0.35 

0.4055 

0.44906 

Long.  Stratified 

0.50909 

0.34927 

0.46119 

Corner  Stratified 

0.57157 

0.43958 

0.40031 

Gaussian 

0.55932 

0.52965 

0.5919 

Gaussian  Side 

0.71881 

0.55738 

0.65176 

Gaussian  Corner 

0.79975 

0.63655 

0.54999 

Total  probability  of  detection  after  50  glimpses  based  on  various  combinations 
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CHAPTER  3: 
Bayesian  Search 


3.1  Theory 

The  following  sections  outline  the  basic  concepts  and  theory  behind  the  Bayesian  approach 
discrete  search  theory. 

3.1.1  Terms 

The  same  terms  from  Chapter  2  apply  along  with  the  following  random  variables: 

Li  -  Target  is  located  in  box  i 
Mi  -  Target  is  not  detected  in  box  i 

po  -  P[tgt  not  in  A]  probability  of  target  not  in  search  area 

n 

By  the  concept  of  total  probability,  £  Pi  -  Pa  and  1  -  Pa  -  Po- 

i= 1 


3.1.2  General  Concepts 

We  use  the  same  basic  assumptions  of  Section  2.1.2  with  an  added  assumption  that  the 
target  may  not  be  in  the  search  area.  The  basics  of  Bayesian  search  theory  is  to  update  the 
probability  of  the  target  being  located  in  the  /th  (p, )  box  based  on  previous  failures  to  detect. 
Searching  the  ith  box  and  failing  to  find  it  will  decrease  its  likelihood  of  being  located  there 
and  increase  the  probability  of  its  being  located  in  another  box.  This  posterior  probability 
comes  from  Bayes’  Theorem. 


P(L,  |  Mi)  = 


P(Mj  1  Lj)P(Lj) 
P{Mi) 


(3.1) 


The  probability  of  not  detecting  the  target  in  the  ith  box  (P(Mj))  is  currently  the  only 
uncalculated  quantity.  It  is,  by  total  probability,  the  sum  of  the  conditional  detection 
failures  based  on  the  location  probabilities.  P(Mj)  =  P(Mj  \  Lj)P(Lj)  +  P(M)  \  -iL,)P(-iL/). 
P(Lj)  and  -i P(Li)  are  pl  and  1  respectively.  P(Mi  \  L;)  =  cp  and  P(Mj  \  -> L,)  =  1  based 
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on  the  assumption  of  now  false  positives.  This  updates  Equation  3.1  to 


Pi  =  P(Lt  |  Mi)  = 


QiPi 

(1  -  pi)  +  pm' 


(3.2) 


The  probability  that  the  target  is  located  in  /th  box  (J  ±  i )  must  also  be  updated  since  it  is 
more  likely  that  it  is  not  located  in  the  /th  box, 


P(Lj  |  Mi) 


P(Mj  1  Lj)P(Lj) 
P(Mi) 


(3.3) 


where  P(L,)  -  pj  and  P{Mf)  =  (1  -  Pi)  +  PM  as  found  above.  P(M)  \  Lf)  is  the  conditional 
probability  of  failing  to  detect  the  target  in  box  /  given  it  is  in  box  /.  Since  there  are  no  false 
positives,  the  searcher  will  always  fail  to  detect  when  searching  in  a  box  where  the  target  is 
not  located,  making  P( M,  \  Lf)  =  1.  This  updates  the  location  probability  by 


P(Lj  I  Mf 


Pj 


(1  -  Pi)  +  PM 


(3.4) 


This  updates  all  pi+j  including  pq. 


3.1.3  Model  Overview 

In  the  Bayesian  search  model,  the  conditional  probabilities  of  detection  failure  only  factor 
into  updating  the  new  location  probabilities.  As  such,  the  search  algorithm  searches  the 
box  with  the  highest  likelihood  of  containing  the  target  and  searches.  Upon  failure  to  find, 
it  updates  all  location  probabilities  and  repeats  until  the  target  is  found  or  the  number  of 
glimpses  has  been  reached. 

This  still  follows  the  Bayesian  search  algorithm  from  Chapter  2  in  that  at  each  step  the 
probability  to  detect  Poi  =  pf  1  -  qf)  is  calculated  for  all  boxes.  The  search  is  conducted  in 
the  box  with  the  highest  detection  probability.  After  a  failure  to  detect,  the  pf  s  are  updated 
via  Equation  3.2  and  Equation  3.4  [14].  The  Bayesian  search  algorithm  then  treats  the  next 
glimpse  as  the  first  with  the  updated  location  probabilities  [14]. 

If  pq  =  0,  this  will  provide  the  same  output  as  the  Discrete  Search  Method.  For  example, 
the  updated  location  probabilities  are  just  A,  from  the  Discrete  Method  divided  by  1  -  <r/,\s 
and  then  normalized  so  that  the  sum  is  equal  to  1.  Multiplying  1  -  qi  s  back  in  after 
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normalization  will  make  the  same  evaluation  as  in  Table  2.1  [14]. 


3.1.4  Examples  from  Other  Research 

The  Bayesian  method  has  been  used  for  several  real-world  searches.  Famously  used  in  USS 
Scorpion,  a  submarine  lost  in  1968,  it  was  also  used  in  MV  Derbyshire,  the  SS  Central 
America,  the  lost  hydrogen  bomb  from  a  B-52  crash.  Air  France  flight  447,  and  most 
recently  flight  MH370.  Figure  3.1  is  the  location  probability  density  originally  constructed 
for  the  USS  Scorpion  [15]. 


Figure  3.1.  USS  Scorpion  Location  Probability.  Source:  [15]. 
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The  basic  steps  outlined  in  previous  examples  are: 
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1.  Construct  plausible  conjectures  as  to  the  cause  of  the  loss  of  the  object. 

2.  For  each  conjecture,  construct  a  location  probability  density  map  (P(L;)  =  /?,). 

3.  Build  the  conditional  detection  probability  maps  for  every  location  covered  under  the 
location  probability  map  ( P(Mj  \  Li)  =  qi).  This  will  be  a  function  of  sensor’s  ability 
to  detect  in  the  given  environment  and  must  be  built  individually  based  on  current 
conditions  and  sensor  dynamics. 

4.  Combine  the  maps  in  Step  2  and  3  to  create  an  overall  detection  probability  density 
map  0,(1  -  qi)). 

5.  Develop  a  search  path  that  scans  the  highest  to  lowest  probability  areas. 

6.  Revise  the  location  probabilities  continuously  during  the  search.  If  multiple  conjec¬ 
tures  are  being  tested  simultaneously,  each  should  be  updated  independently  and  if 
the  search  fails  to  corroborate  aspects  of  a  particular  conjecture  the  likelihood  of  that 
scenario  will  need  to  be  downgraded  [15]. 

One  last  aspect  of  the  Bayesian  search,  that  is  not  being  explored  in  this  research,  is  the 
concept  of  the  cost  of  a  search.  If,  for  example,  searching  one  particular  box  takes  longer 
time,  requires  more  resources,  or  is  more  dangerous,  then  that  should  be  taken  into  account 
as  a  weight  on  the  search.  If  the  cost  function  can  be  estimated  as  q,  then  the  new  detection 
probabilities  come  out  as  Pl^~q,)  [16], 

3.2  Problem  and  Model 

This  basic  model  for  Bayesian  search  is  based  solely  on  the  conditional  miss  probabilities 
(qi)  and  the  initial  location  probabilities  (p0i ).  These  parameters  will  be  varied  and  results 
shown. 

3.2.1  Parameters 

The  location  and  conditional  miss  will  have  the  same  variations  as  presented  in  Chapter  2 
(Figures  2.1  and  2.2)  and  they  are  repeated  in  Figures  3 .2  and  3 .3  for  the  reader’s  convenience 
with  po  =  .1  and  the  detection  probabilities  will  be  varied  according  to  Figure  2.2,  which 
are  shown  again  here  for  clarity. 
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Figure  3.2.  Bayesian  Initial  Location  Probabilities  p0i 


From  top  left  clockwise:  Uniform  dist.,  Linear  stratified  along  longitude,  Linear 
stratified  along  Long/Lat,  Gaussian  dist.,  Gaussian  centered  in  corner,  Gaussian 
centered  on  side 


Figure  3.3.  Bayesian  Miss  Probabilities 
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3.3  Results 

The  following  sections  show  the  results  for  six  different  location  probability  densities  and 
the  three  conditional  search  probability  densities.  Each  combination  was  examined  for 
n  =  50  glimpses  (half  the  number  of  available  boxes). 

The  tables  and  figures  showing  all  18  possible  combinations  are  presented  in  Appendix  A 
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and  Appendix  B,  respectively. 


3.3.1  Linear  Stratified  Location 

Figures  3.4,  3.5,  and  3.6  show  the  remaining  location  probability  densities  for  the  three 
conditional  search  probability  densities  after  50  glimpses.  Also  shown  are  the  boxes  that 
were  selected  to  search  by  the  Bayesian  search  algorithm. 


For  the  uniform  conditional  search  probability  in  Figure  3.4,  unsurprisingly,  the  selected 
search  boxes  start  with  highest  p,  and  slowly  work  to  the  left  as  p/  lowers.  This  follows 
the  same  pattern  as  the  discrete  search  from  Chapter  2.  The  exception  being  that  the 
upper-leftmost  graph  of  Figure  3.4  now  shows  the  updated  location  probabilities  in  Bayes’ 
algorithm  instead  of  merely  the  next  search  step  as  in  the  discrete  search  method.  This 
follows  as  the  Bayesian  method  and  Discrete  method  calculate  roughly  the  same  probability 
of  detection  with  the  exception  that  the  Bayesian  search  modifies  po  =  .1  — »  0.1846  . 


Figure  3.4.  Stratified  Location  and  Uniform  Search  Density  Results 


Remaining  probability  of  detection  and  boxes  search  based  on  stratified  p,-  and 
uniform 


For  the  longitudinal  stratified  conditional  search  probability  in  Figure  3.5,  the  selected 
search  boxes  also  start  with  highest  p,  and  slowly  work  to  the  left  as  p,-  lowers.  This  follows 
the  slightly  unexpected  results  from  the  discrete  search  of  Chapter  2  possibly  from  the  factor 
of  1  -  qi  being  93.7%  in  the  western  side  and  37%  in  the  east.  This  coupled  with  the  p; 
being  a  factor  of  10  higher  in  the  east  than  the  west  leads  it  to  be  a  higher  probability  of 
detection  even  though  the  sensor  performs  worse  there.  The  probability  of  the  target  not 
being  in  the  search  area  p0  =  .1  — >  0.1458. 
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Figure  3.5.  Stratified  Location  and  Longitudinal  Stratified  Search  Results 
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For  the  corner  stratified  conditional  search  probability  in  Figure  3.6,  the  selected  boxes 
start  in  the  southeast  due  to  the  juxtaposition  of  high  pt  and  low  qt.  The  path  moves 
mostly  along  discretized  radial  arcs  centered  on  the  southeast  corner.  Which  box  is  chosen 
via  the  Bayesian  search  algorithm  is  dependent  on  when  1  -  qt  will  push  a  probability  of 
detection  higher  for  a  lower  pj.  The  probability  of  the  target  not  being  in  the  search  area 
po  =  .1  ->  0.1710. 


Figure  3.6.  Stratified  Location  and  Corner  Stratified  Search  Density  Results 


Remaining  probability  of  detection  and  boxes  search  based  on  stratified  p,  and 
corner  stratified  q{ 


For  the  longitudinal  stratified  p,  location  probability  density,  the  following  results  were 
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found  based  on  running  the  search  for  up  to  500  glimpses  or  a  total  probability  of  detection 
at  95%.  A  data  set  using  a  step  size  of  5%  is  plotted  in  Figure  3.7  along  with  associated  fit 
lines  as  calculated  from  MATLAB’s  fit  algorithm  based  on 

fix)  =  a  -  e~bx. 


Each  of  the  fit  lines  is  color  coded  to  the  data  set  with  black  for  uniform,  red  for  longitudinal 
stratification,  and  green  for  corner  stratification.  As  the  number  of  glimpses  increases,  the 
probability  approaches  Pa  =  .9,  with  the  uniform  case  increasing  faster  than  the  other  two. 
The  associated  parameters  for  the  best  fit  lines  are  presented  in  Table  3.1. 

Figure  3.7.  Stratified  Location  Cumulative  Probability  and  Fit  Lines 


F*d  vs  Number  of  Glimpses  and  Exponential  Fitted  Curves 


Cumulative  probabilities  of  detection  based  on  the  number  of  glimpses.  Included 
are  the  fit  best  fit  lines 


The  results  are  very  similar  to  results  from  Chapter  2  as  the  same  decision  algorithm  is 
used;  however,  the  exponential  growth  approaches  Pa  =  -9  the  initial  probability  that  the 
target  is  in  the  search  area  A. 
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Table  3.1.  Parameter  Variation  for  p,  for  Longitudinal  Stratification 


Parameter 

a 

b 

r2 

SSE 

RMSE 

Pstrat  ^uniform 
Pstratitflcmg 
Pstrat  incomer 

0.8938 

0.90437 

0.89101 

-0.015155 

-0.0097729 

-0.0131 

0.99013 

0.99496 

0.97996 

0.0086678 

0.0053481 

0.017876 

0.0093465 

0.017237 

0.031514 

Best  fit  line  parameters  for  a  longitudinally  stratified  location  probability  density 
of  the  form  f(x)  =  a  -  e~bx 


3.3.2  Gaussian  Corner  Location 

Figures  3.8,  3.9,  and  3.10  show  the  remaining  location  probability  densities  for  the  three 
conditional  search  probability  densities  after  50  glimpses.  Also  shown  are  the  boxes  that 
were  selected  to  search  by  the  Bayesian  search  algorithm. 


For  the  uniform  conditional  search  probability  in  Figure  3.8,  unsurprisingly,  the  selected 
boxes  started  with  highest  pi  and  slowly  works  outward  as  pi  shifts  higher  probability  to 
other  boxes.  A  second  look  begins  occurring  in  those  higher  probability  boxes  prior  to 
searching  all  boxes. 


Figure  3.8.  Gaussian  Corner  Location  and  Uniform  Search  Density  Results 


Remaining  probability  of  detection  and  boxes  search  based  on  Gaussian  pi 
and  uniform  qi 


For  the  longitudinal  stratified  conditional  search  probability  in  Figure  3.9,  the  selected  boxes 
also  started  with  highest  pt  and  slowly  work  outward  as  pi  lowers.  This  follows  the  slightly 
unexpected  result  from  Chapter  2  since  there  are  better  probabilities  of  detection  in  some 
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areas  of  lower  target  location  probability.  As  the  location  probability  shifts  to  other  boxes, 
the  higher  pi  increases  the  likelihood  of  the  target  being  located  to  the  west  of  the  initial 
location  locus.  From  the  Remainder  Probability  graph  in  Figure  3.9,  the  higher  location 
probability  remains  dispersed  in  the  northwest  corner  albeit  non- uniformly. 


Figure  3.9.  Gaussian  Corner  and  Longitudinal  Stratified  Search  Density  Re¬ 
sults 


Remaining  probability  of  detection  and  boxes  search  based  on  stratified  p,  and 
linear  stratified  qi 


For  the  corner  stratified  conditional  search  probability  in  Figure  3.10,  the  selected  boxes 
start  in  the  northeast  due  to  the  high  pi  overriding  the  high  q,.  The  Remainder  Probability 
graph  in  Figure3.10  shows  some  areas  2-3  times  greater  than  other  areas,  especially  around 
the  original  locus.  Some  boxes  have  been  searched  twice  or  three  times  before  every  box 
has  been  searched  once.  Since  the  location  probability  distribution  pi  remains  identical  to 
that  from  Chapter  2  (other  than  the  being  90%  of  the  original),  the  overall  patterns  of  the 
search  are  identical  to  the  results  from  Chapter  2. 
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Figure  3.10.  Gaussian  Corner  and  Corner  Stratified  Search  Density  Results 


Remaining  probability  of  detection  and  boxes  search  based  on  stratified  p,-  and 
corner  stratified  g,; 


For  the  Gaussian  corner  pi  location  probability  density,  the  following  results  were  found 
based  on  running  the  search  for  up  to  500  glimpses  or  a  total  probability  of  detection  at 
95%.  A  data  set  using  a  step  size  of  5%  is  plotted  in  Figure  3.11  along  with  associated  fit 
lines  as  calculated  from  MATLAB’s  fit  algorithm  based  on 

f(x)  =  a  -  e~b\ 


Each  of  the  fit  lines  are  color  coded  to  the  data  set  with  black  for  uniform,  red  for  longitudinal 
stratification,  and  green  for  corner  stratification.  As  the  number  of  glimpses  increase,  the 
probability  approaches  Pa  =  .9,  with  the  uniform  case  increasing  faster  than  the  other  two. 
The  associated  parameters  for  the  best  fit  lines  are  presented  in  Table  3.2. 
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Figure  3.11.  Gaussian  Corner  Cumulative  Probability  and  Fit  Lines 


F*d  vs  Number  of  Glimpses  and  Exponential  Fitted  Curves 


Cumulative  probabilities  of  detection  based  on  the  number  of  glimpses.  Included 


are  the  fit  best  fit  lines 


Table  3.2.  Parameter  Variation  for  p,  for  Corner  Gaussian  Distribution 


Parameter 

a 

b 

71 

SSE 

RMSE 

PGC  iQ uniform 
PGCiQlong 
PGCi<1  corner 

0.89231 

0.88774 

0.88725 

-0.03706 

-0.022588 

-0.01732 

0.98424 

0.97697 

0.9747 

0.0068884 

0.014581 

0.018924 

0.019563 

0.028462 

0.032424 

Best  fit  line  parameters  for  a  corner  Gaussian  location  probability  density  of  the 
form  fix)  =  a  -  e~bx 


The  results  are  very  similar  to  results  from  Chapter  2  as  the  same  decision  algorithm  is 
used. 


3.3.3  Cumulative  Detection  Probability 

The  total  probability  of  detection  produced  by  the  Bayesian  search  algorithm  after  50 
glimpses  is  shown  in  Table  3.3.  The  Gaussian  p,-’  s  have  higher  probabilities,  which  is  not 
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unexpected  since  they  have  a  much  higher  concentration  on  a  few  points.  This  allows  for 
multiple  searches  over  a  smaller  number  of  points  instead  of  50  different  boxes.  The  other 
result  seems  to  be  that  higher  probability  results  when  p,  matches  up  to  1  -  cp  as  it  does 
when  p stratified  has  a  higher  probability  of  detection  when  paired  with  qComerstrat  instead 
of  qiongstrat-  This  follows  the  same  pattern  as  the  Discrete  Search  Method  but  lower  overall 
since  there  is  some  probability  that  the  target  is  outside  of  A.  The  increase  to  po  after  50 
glimpses  is  shown  in  Table  3.4. 


Table  3.3.  Pd  after  50  glimpses 


Pi 

Uniform  Long.  Stratified 

Corner  Stratified 

Uniform 

0.315 

0.36495 

0.40415 

Long.  Stratified 

0.45818 

0.31435 

0.41507 

Corner  Stratified 

0.51441 

0.39562 

0.36028 

Gaussian 

0.50339 

0.47669 

0.53271 

Gaussian  Side 

0.64693 

0.50164 

0.58658 

Gaussian  Corner 

0.71978 

0.57289 

0.49499 

Total  probability  of  detection  after  50  glimpses  based  on  various  combinations 


Table  3.4.  po  after  50  glimpses 


x\  qt 

Pi  \ 

Uniform  Long.  Stratified 

Corner  Stratified 

Uniform 

0.1460 

0.1575 

0.1678 

Long.  Stratified 

0.1846 

0.1458 

0.1710 

Corner  Stratified 

0.2059 

0.1655 

0.1563 

Gaussian 

0.2014 

0.1911 

0.2140 

Gaussian  Side 

0.2832 

0.2007 

0.2419 

Gaussian  Corner 

0.3569 

0.2341 

0.1980 

Probability  of  target  not  in  A  after  50  glimpses  based  on  various  combinations 
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CHAPTER  4: 
Continuous  Search 


4.1  Theory 

The  following  section  outlines  the  basic  concepts  and  theory  behind  a  continuous  exhaustive 
search. 

Terms 

A  -  Total  search  area  (no  searches  conducted  outside  of  A) 
vs  -  Searcher’s  velocity 
vt  -  Target’s  velocity 

R  -  Range  of  the  sensor.  Cookie  cutter  sensor  that  will  detect  the  target  if  the  range 
is  less  than  R.  R  =  R(vs) 
w  -  Sweep  width,  w  =  2R 

The  sensor  is  a  "cookie  cutter"  sensor  that  has  100%  chance  of  detecting  the  target  if  the 
distance  from  the  searcher  to  the  target  is  less  than  R.  In  [6]  both  R  and  both  vs  were 
independent  of  each  other  and  varied  accordingly.  Here,  R  is  a  function  of  vs  which  more 
closely  resembles  real-world  sensor  performance,  especially  in  SONAR  sensors. 

4.1.1  General  Concepts 

The  basic  assumptions  for  this  continuous  search  model  are  [6]  [9]  [13]: 

•  Single  target  undergoing  Brownian  diffusion 

•  Two  searchers  with  "cookie-cutter"  sensors 

•  Exhaustive  search  using  one  of  three  search  patterns 

•  Sensors  radius  is  velocity  dependent  R(vs ) 

•  Each  searcher  results  independent  of  the  other 

•  No  false  positives 

•  Non-adaptive  search. 
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As  shown  in  Figure  4.1  the  sweep  width  for  a  moving  searcher  (w  =  2R )  will  cover  an 
area  searched  A{t )  =  wvt.  The  exhaustive  search  covers  all  of  A  without  overhanging 
any  previously  searched  area  or  any  area  outside  of  A.  The  total  time  to  coverage  then  is 


Figure  4.1.  Search  Area  and  Sweep  Width.  Source:  [9]. 


Search  Area  (A)  as  it  relates  to  Sweep  Width  (u>) 


4.1.2  Model  Overview 

Based  on  the  work  done  by  H.  Wang  and  H.  Zhou  [9]  and  Z.  Lukens’s  thesis  [6],  the  two 
searchers  will  move  along  one  of  three  set  paths:  the  inline  ladder,  the  multipath  ladder  and 
the  inline  spiral  search.  For  each  path  setup,  100,000  independent  Monte  Carlo  simulations 
are  conducted  with  the  target  uniformly  distributed  over  A.  The  independent  targets  move 
randomly  with  a  Brownian  motion  with  a  given  diffusion  rate. 

The  searchers  will  move  along  the  given  path  independently  of  each  other  without  coordi¬ 
nation.  An  input  parameter  of  the  initial  separation  of  the  two  searcher  is  given  so  that  the 
searchers  don’t  overlap  throughout  the  search.  The  separation  in  the  inline  searches  was 
estimated  at  4LX  where  Lx  is  half  the  width  of  A  [6]. 

4.2  Problem  and  Model 

The  model  is  thus  two  searchers  (aircraft,  surface  vessels,  subs,  etc.)  of  similar  capability 
performing  the  inline  spiral  in  Figure  4.2,  the  inline  ladder  search  in  Figure  4.3,  and  the 
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multipath  ladder  search  in  Figure  4.4. 


Figure  4.2.  Spiral  Search.  Source:  [9]. 


Figure  4.3.  Inline  Ladder  Search.  Source:  [9]. 


Figure  4.4.  Multipath  Ladder  Search  Pattern.  Source:  [9]. 


4.2.1  Parameter 

As  an  extension  of  previous  work  done  [6],  this  chapter  will  explore  the  search  patterns 
with  a  sensor  that  obeys  Equation  4.1. 

R(VS)  =  R0e~Vs/v  (4.1) 

Rq  is  the  performance  of  the  sensor  at  vs  =  0  and  V  is  the  exponential  decay  factor  based 
on  the  sensor  performance  at  speed.  Each  sensor,  in  reality,  will  behave  differently  based 
on  its  characteristics.  For  instance,  a  towed  SONAR  array  requires  a  minimum  speed  in 
order  to  be  drooped  correctly  and  at  the  proper  depth,  whereas  a  hull-mounted  SONAR 
does  not  have  this  minimum  effect  and  is  only  limited  at  higher  speed  by  self-noise.  Further 
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compounding  the  real  life  sensors  is  that,  at  higher  speeds,  the  self  noise  has  non-linear 
effects  based  on  cavitation  above  a  threshold.  To  simplify  the  issue,  the  sensor  is  assumed 
to  be  operating  in  the  range  where  Equation  4.1  holds.  Each  search  velocity  will  be  varied 
as  vs  6  (8, 9, 10, 11, 12)  and  the  velocity  decay  factor  from  V  e  (5,  6, 7,  8, 9, 10)  giving  150 
different  combinations  for  each  search  pattern.  The  initial  search  radius,  Rq  =  2.5  nm,  will 
be  set  based  on  the  largest  search  radius  in  [6]. 


4.3  Results 

In  Tables  4.1,  4.2,  and  4.3,  the  time  in  hours  to  achieve  a  probability  of  detection  at  90% 
(Pf)  =  90%)  are  shown  for  V  =  9  the  inline  ladder  search,  the  multipath  ladder,  and  the 
spiral  search,  respectively. 

Table  4.1.  Time  to  Po  =  90%  (V  =  9) 


Velocities 

8 

9 

Ship  1 

10 

11 

12 

8 

6.8662 

6.9782 

7.1117 

7.5250 

7.7083 

9 

7.5069 

7.8241 

7.7043 

7.8639 

7.9255 

Ship  2  10 

9.7362 

9.0589 

8.7202 

8.4741 

8.3054 

11 

10.1011 

10.0333 

9.2198 

8.9486 

8.6854 

12 

10.0425 

10.1511 

10.0968 

9.2826 

8.9569 

Inline  ladder  search 
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Table  4.2.  Time  to  Pd  =  90%(V  -  9) 


Velocities 

8 

9 

Ship  1 

10 

11 

12 

8 

8.1866 

8.2470 

8.2123 

8.1351 

8.0883 

9 

9.3043 

8.9871 

8.7202 

8.3385 

8.1969 

Ship  2  10 

10.3288 

9.5669 

9.2282 

8.7452 

8.4140 

11 

10.7790 

10.0333 

9.4231 

9.0164 

8.6854 

12 

10.4768 

10.5311 

9.7168 

9.1197 

8.7940 

Multipath  ladder  search 


Table  4.3.  Time  to  Pd  =  90%(V  =  9) 


Velocities 

8 

9 

Ship  1 

10 

11 

12 

8 

9.6390 

9.3043 

9.6515 

10.2367 

10.3139 

9 

9.8330 

9.4100 

9.3129 

9.7621 

9.8254 

Ship  2  10 

10.4135 

9.4822 

9.3129 

9.5587 

9.6625 

11 

10.5078 

10.0333 

9.5587 

9.6943 

9.7168 

12 

10.3682 

11.6710 

11.5625 

10.6397 

10.5854 

Spiral  search 


In  a  marked  difference  from  the  results  found  in  [6],  faster  search  speed  does  not  result  in  a 
quicker  time  to  reach  Pd  =  ■ 9  in  all  cases.  For  instance,  the  lowest  time  to  detection  for  the 
inline  ladder  search  is  with  v\  =  8  and  vo  =  8,  for  multipath  ladder  v\  =  12  and  vo  =  8,  and 
for  spiral  search  v\  =  9  and  V2  =  8.  In  fact,  keeping  vo  =  8  and  varying  v\  over  the  range 
causes  an  increase  and  then  decrease.  Similar  results  happen  on  some  other  variations  but 
as  can  be  seen  it  is  not  consistent  from  one  speed  to  another.  Another  noteworthy  result  is 
that  varying  the  speed  of  ship  1  has  less  of  an  effect  normally  than  varying  ship  2’s  speed. 
An  example  of  this  can  be  seen  in  the  multipath  ladder.  Raising  v\  from  8  to  12  knots  keeps 
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the  time  to  detection  relatively  constant  around  8  hours,  while  raising  v\  from  8  to  12  knots 
increase  the  time  to  detect  to  approximately  8.5  hours. 

Similar  to  the  results  found  in  [6],  the  inline  ladder  search  has  the  lowest  time  to  detect, 
followed  by  the  multipath  ladder,  and  finally  the  spiral  search  for  a  given  set  of  parameters. 
As  expected,  increasing  V  only  lowers  the  time  to  detection.  This  follows  [6]  as  increasing 
V  raises  the  search  radii  without  changing  the  speed.  Direct  comparison  to  other  results 
from  [6]  is  problematic.  This  work  varies  three  of  the  four  parameters  from  Equation  4.1, 
directly  varying  v\  and  vo  and  indirectly  R\  and  Ri.  In  [6],  R\  and  Ro  are  varied  directly.  If 
Equation  4.1  is  constructed  with  a  V  and  Rq  so  that  R\  and  /sS  equal  values  from  [6],  then 
results  will  be  the  same  as  the  search  algorithm  work  in  the  same  manner. 

A  similar  presentation  for  the  various  searches  is  outlined  below  in  Figures  4.5-4.10. 

Spiral  Search 

Figure  4.5.  Pd  vs  Time  Figure  4.6.  Pd  vs  Time 
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Multipath  Ladder 

Figure  4.7.  Pd  vs  Time  Figure  4.8.  Pd  vs  Time 


0  5  10  15  20  25  0  5  10  15  20  25 

Time  (hrs)  Time  (hrs) 


Inline  Ladder  Search 


Figure  4.9.  Pd  vs  Time 

Figure  4.10.  Pd  vs  Time 

R, =0.599,  R2=0.599,  V(=10,  V2=10,  V=7 
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As  can  be  seen  from  the  figures,  there  is  an  approach  to  approximately  Pd  -  1  -  e~bt .  An 
increase  in  V,  as  stated  before  increases  R(vs)  for  a  given  vs.  This  obviously  hastens  PD  to 
higher  level  more  quickly. 
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CHAPTER  5: 

Conclusions  and  Future  Work 


Based  on  the  results  from  Chapters  2-4  the  following  sections  summarize  the  conclusions  drawn. 

5.1  Results  and  Conclusions 

5.1.1  Discrete  and  Bayesian  Searches 

Both  the  Discrete  and  Bayesian  searches  follow  the  greedy  algorithm  and  give  similar 
results.  In  fact,  if  po  =  0,  then  the  results  will  be  identical  to  within  rounding  error.  As 
expected,  the  more  concentrated  p,-  density,  the  higher  the  probability  to  detect  Pd  and  the 
lower  number  of  glimpse  needed  to  reach  it.  This  even  has  the  capacity  to  override  higher 
conditional  miss  probability  in  that  same  concentrated  area.  Both  search  systems  provide 
predictable,  repeatable  results,  but  the  Bayesian  search  provides  a  more  useful  output.  By 
modifying  the  p,  at  each  step,  not  only  can  po  >  0,  but  faulty  conjectures  that  created  the 
Pi  be  identified  by  scrutinizing  the  changing  pi  and  increasing  po-  This  provides  a  useful 
product  by  the  algorithm  at  any  point,  the  p,  have  been  calculated  if  the  number  of  glimpses 
have  been  exceeded  and  can  be  used  at  a  later  date.  For  a  discrete  search,  the  output  from 
the  algorithm  is  just  the  failure  to  detect  on  previous  searches. 

5.1.2  Continuous  Searches 

As  vs  increases,  Pd  has  a  non-intuitive  result.  It  varies  in  a  non-linear,  way  some  times 
increasing  or  decreasing  as  vs  changes.  Another  non-intuitive  result  is  that  v\  has  a  different 
change  in  PD  than  vo.  The  last  result,  an  expected  one,  is  that  increasing  V  increases  the  PD 
by  increasing  R. 


5.2  Future  Work 

5.2.1  Discrete  and  Bayesian  Searches 

Improvements  to  the  Discrete  and  Bayesian  model  in  future  work  could  include: 
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•  Multiple  pi  tested  simultaneously  updating  all  based  on  multiple  conjectures 

•  Multiple  <27  tested  simultaneously  updating  all  based  on  multiple  exclusive  sensors 

•  Time  dependent  such  as  c/j  =  e~r,/T  as  greedy  algorithm  not  the  optimal  solution 

•  Time  dependent  pi  for  a  moving  target. 


5.2.2  Continuous  Searches 

Improvements  to  the  continuous  search  in  future  work  could  include: 

•  Different  search  patterns 

•  Non-cookie  cutter  sensor  with  some  non-uniform  probability  of  detection 

•  Game  theory  approach  with  a  target  counter-detecting  and  evading. 
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APPENDIX  A: 
Figures 


A.l  Discrete 


Uniform  Location  and  Uniform  Search  Density  Results 


Figure  A.l.  Remaining  probability  of  detection  and  boxes  search  based  on 
uniform  /?,■  and  uniform  qi 


Uniform  Location  and  Longitudinal  Stratified  Search  Density  Results 


Figure  A. 2.  Remaining  probability  of  detection  and  boxes  search  based  on 
uniform  /?,  and  linear  stratified  <37 
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Uniform  Location  and  Corner  Stratified  Search  Density  Results 


Figure  A. 3.  Remaining  probability  of  detection  and  boxes  search  based  on 
uniform  p\  and  corner  stratified  q\ 


Uniform  Location  Cumulative  Probability  and  Fit  Lines 


F*d  vs  Number  of  Glimpses  and  Exponential  Pitted  Curves 


Figure  A. 4.  Cumulative  probabilities  of  detection  based  on  the  number  of 
glimpses.  Included  are  the  fit  best  fit  lines 
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Corner  Stratified  Location  and  Uniform  Search  Density  Results 


Figure  A. 5.  Remaining  probability  of  detection  and  boxes  search  based  on 
corner  stratified  p,  and  uniform  qi 


Corner  Stratified  Location  and  Longitudinal  Stratified  Search  Density  Results 


Figure  A. 6.  Remaining  probability  of  detection  and  boxes  search  based  on 
corner  stratified  p;  and  linear  stratified 
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Corner  Stratified  Location  and  Corner  Stratified  Search  Density  Results 


Figure  A. 7.  Remaining  probability  of  detection  and  boxes  search  based  on 
corner  stratified  pi  and  corner  stratified  qi 


Corner  Stratified  Location  Cumulative  Probability  and  Fit  Lines 


F*d  vs  Number  of  Glimpses  and  Exponential  Pitted  Curves 


Figure  A. 8.  Cumulative  probabilities  of  detection  based  on  the  number  of 
glimpses.  Included  are  the  fit  best  fit  lines 
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Gaussian  Location  and  Uniform  Search  Density  Results 


Figure  A. 9.  Remaining  probability  of  detection  and  boxes  search  based  on 
Gaussian  and  uniform  q\ 


Gaussian  Location  and  Longitudinal  Stratified  Search  Density  Results 


Figure  A. 10.  Remaining  probability  of  detection  and  boxes  search  based  on 
Gaussian  p ,■  and  linear  stratified  qi 
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Gaussian  Location  and  Corner  Stratified  Search  Density  Results 


Figure  A. 11.  Remaining  probability  of  detection  and  boxes  search  based  on 
Gaussian  and  corner  stratified  q\ 


Gaussian  Location  Cumulative  Probability  and  Fit  Lines 


Figure  A. 12.  Cumulative  probabilities  of  detection  based  on  the  number  of 
glimpses.  Included  are  the  fit  best  fit  lines 
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Side  Gaussian  Location  and  Uniform  Search  Density  Results 


Figure  A. 13.  Remaining  probability  of  detection  and  boxes  search  based  on 
side  Gaussian  p ,■  and  uniform  qi 


Side  Gaussian  Location  and  Longitudinal  Stratified  Search  Density  Results 


Long. 


Long. 


Figure  A. 14.  Remaining  probability  of  detection  and  boxes  search  based  on 
side  Gaussian  p,-  and  linear  stratified  qt 
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Side  Gaussian  Location  and  Corner  Stratified  Search  Density  Results 


Figure  A. 15.  Remaining  probability  of  detection  and  boxes  search  based  on 
side  Gaussian  p ,■  and  corner  stratified  q\ 


Side  Gaussian  Location  Cumulative  Probability  and  Fit  Lines 


Figure  A. 16.  Cumulative  probabilities  of  detection  based  on  the  number  of 
glimpses.  Included  are  the  fit  best  fit  lines 
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A.2  Bayesian 

Uniform  Location  and  Uniform  Search  Density  Results 


Figure  A. 17.  Remaining  probability  of  detection  and  boxes  search  based  on 
uniform  pi  and  uniform  q\ 


Uniform  Location  and  Longitudinal  Stratified  Search  Density  Results 


Figure  A. 18.  Remaining  probability  of  detection  and  boxes  search  based  on 
uniform  /?,■  and  linear  stratified  <7/ 
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Uniform  Location  and  Corner  Stratified  Search  Density  Results 


Figure  A. 19.  Remaining  probability  of  detection  and  boxes  search  based  on 
uniform  /?,■  and  corner  stratified  q, ■ 


Uniform  Location  Cumulative  Probability  and  Fit  Lines 


F*d  vs  Number  of  Glimpses  and  Exponential  Pitted  Curves 


Figure  A. 20.  Cumulative  probabilities  of  detection  based  on  the  number  of 
glimpses.  Included  are  the  fit  best  fit  lines 
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Corner  Stratified  Location  and  Uniform  Search  Density  Results 


Figure  A. 21.  Remaining  probability  of  detection  and  boxes  search  based  on 
corner  stratified  p /  and  uniform  qj 


Corner  Stratified  Location  and  Longitudinal  Stratified  Search  Density  Results 


Figure  A. 22.  Remaining  probability  of  detection  and  boxes  search  based  on 
corner  stratified  pi  and  linear  stratified 
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Corner  Stratified  Location  and  Corner  Stratified  Search  Density  Results 


Figure  A. 23.  Remaining  probability  of  detection  and  boxes  search  based  on 
corner  stratified  p ;  and  corner  stratified  qi 


Corner  Stratified  Location  Cumulative  Probability  and  Fit  Lines 


F*d  vs  Number  of  Glimpses  and  Exponential  Pitted  Curves 


Figure  A. 24.  Cumulative  probabilities  of  detection  based  on  the  number  of 
glimpses.  Included  are  the  fit  best  fit  lines 


52 


Gaussian  Location  and  Uniform  Search  Density  Results 


Figure  A. 25.  Remaining  probability  of  detection  and  boxes  search  based  on 
Gaussian  and  uniform  q\ 


Gaussian  Location  and  Longitudinal  Stratified  Search  Density  Results 


Figure  A. 26.  Remaining  probability  of  detection  and  boxes  search  based  on 
Gaussian  p ,■  and  linear  stratified  qi 
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Gaussian  Location  and  Corner  Stratified  Search  Density  Results 


Figure  A. 27.  Remaining  probability  of  detection  and  boxes  search  based  on 
Gaussian  and  corner  stratified  q\ 


Gaussian  Location  Cumulative  Probability  and  Fit  Lines 


Figure  A. 28.  Cumulative  probabilities  of  detection  based  on  the  number  of 
glimpses.  Included  are  the  fit  best  fit  lines 
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Side  Gaussian  Location  and  Uniform  Search  Density  Results 


Figure  A. 29.  Remaining  probability  of  detection  and  boxes  search  based  on 
side  Gaussian  p ,■  and  uniform  qi 


Side  Gaussian  Location  and  Longitudinal  Stratified  Search  Density  Results 


Figure  A. 30.  Remaining  probability  of  detection  and  boxes  search  based  on 
side  Gaussian  p,-  and  linear  stratified  qt 
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Side  Gaussian  Location  and  Corner  Stratified  Search  Density  Results 
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Figure  A. 31.  Remaining  probability  of  detection  and  boxes  search  based  on 
side  Gaussian  p ,■  and  corner  stratified  q\ 


Side  Gaussian  Location  Cumulative  Probability  and  Fit  Lines 


F*d  vs  Number  of  Glimpses  and  Exponential  Pitted  Curves 


Figure  A. 32.  Cumulative  probabilities  of  detection  based  on  the  number  of 
glimpses.  Included  are  the  fit  best  fit  lines 
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Figure  A. 97.  Pd  vs  Time  Figure  A. 98.  Pd  vs  Time 


5  10  15  20  25  0  5  10  15  20 

Time  (hrs)  Time  (hrs) 


Figure  A. 99.  Pd  vs  Time  Figure  A. 100.  Pd  vs  Time 
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Figure  A. 121.  Pd  vs  Time 
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Figure  A. 123.  Pd  vs  Time 
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Figure  A. 337.  Pd  vs  Time  Figure  A. 338.  Pd  vs  Time 
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Figure  A. 340.  Pd  vs  Time 


133 


Probability  of  Detection  Probability  of  Detection 


Time  (hrs) 

Figure  A. 341.  Pq  vs  Time 

R1  =0.797,  R2=0.797,  V  =8,  V2=8,  V=7 

1  - 1 - ' - T - 

0.9  - 
0.8  - 

0.7  -  / 

0.6  -  / 

0.5  -  / 

0.4  - 


o' - 1 - 1 - 1 - 1 - 

0  5  10  15  20 

Time  (hrs) 

Figure  A. 343.  Pd  vs  Time 


Time  (hrs) 

Figure  A. 342.  Pd  vs  Time 


Time  (hrs) 

Figure  A. 344.  Pd  vs  Time 


134 


Probability  of  Detection  Probability  of  Detection 


Figure  A. 345.  Pd  vs  Time 
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Figure  A. 353.  Pd  vs  Time 
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Figure  A. 360.  Pd  vs  Time 
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Figure  A. 361.  Pd  vs  Time  Figure  A. 362.  Pd  vs  Time 
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Figure  A. 364.  Pd  vs  Time 
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Figure  A. 368.  Pd  vs  Time 
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Figure  A. 371.  Pd  vs  Time 
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Figure  A. 372.  Pd  vs  Time 
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Figure  A. 377.  Pd  vs  Time  Figure  A. 378.  Pd  vs  Time 
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Figure  A. 380.  Pd  vs  Time 


143 


Probability  of  Detection  Probability  of  Detection 


Figure  A. 381.  Pd  vs  Time 


Figure  A. 383.  Pd  vs  Time 
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Figure  A. 384.  Pd  vs  Time 
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Figure  A. 387.  Pd  vs  Time 
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Figure  A. 388.  Pd  vs  Time 
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Figure  A. 389.  Pd  vs  Time 


Figure  A. 390.  Pd  vs  Time 


Figure  A. 391.  Pd  vs  Time 


Figure  A. 392.  Pd  vs  Time 
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Figure  A. 393.  Pd  vs  Time 


Figure  A. 395.  Pd  vs  Time 
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Figure  A. 396.  Pd  vs  Time 
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Figure  A. 397.  Pd  vs  Time 
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Figure  A. 404.  Pd  vs  Time 
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Figure  A. 405.  Pd  vs  Time 


Figure  A. 407.  Pd  vs  Time 
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Figure  A. 408.  Pd  vs  Time 
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Figure  A. 413.  Pd  vs  Time 
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Figure  A. 417.  Pd  vs  Time  Figure  A. 418.  Pd  vs  Time 


Figure  A. 419.  Pd  vs  Time 


Figure  A. 420.  Pd  vs  Time 
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Figure  A. 421.  Pd  vs  Time  Figure  A. 422.  Pd  vs  Time 


Figure  A. 423.  Pd  vs  Time 


Figure  A. 424.  Pd  vs  Time 
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Figure  A. 428.  Pd  vs  Time 
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Figure  A. 431.  Pd  vs  Time 
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Figure  A. 432.  Pd  vs  Time 
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Figure  A. 437.  Pd  vs  Time  Figure  A. 438.  Pd  vs  Time 
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Figure  A. 439.  Pd  vs  Time 


Figure  A. 440.  Pd  vs  Time 
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Figure  A. 441.  Pd  vs  Time 


Figure  A. 443.  Pd  vs  Time 


R, =0.920,  R2=0.558,  V  =8,  V  =12,  V= 8 

1  - T - T - —  [ - 

0.9- 
0.8  - 
|  0.7  - 

I  0.6- 

□ 

o  0.5  -  / 

1  0-4  -  / 

£i 

0  / 

n_  0.3  -  / 

0.2-  / 

0.1  -I 

01 - 1 - 1 - 1 - 1 - 

0  5  10  15  20  25 

Time  (hrs) 

Figure  A. 442.  Pd  vs  Time 


Figure  A. 444.  Pd  vs  Time 
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Figure  A. 445.  Pd  vs  Time  Figure  A. 446.  Pd  vs  Time 


Figure  A. 447.  Pd  vs  Time  Figure  A. 448.  Pd  vs  Time 
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Figure  A. 449.  Pd  vs  Time 


Figure  A. 450.  Pd  vs  Time 


Figure  A. 451.  Pd  vs  Time 


Figure  A. 452.  Pd  vs  Time 
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Figure  A. 453.  Pd  vs  Time 


Figure  A. 455.  Pd  vs  Time 
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Figure  A. 456.  Pd  vs  Time 
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Figure  A. 457.  Pd  vs  Time  Figure  A. 458.  Pd  vs  Time 


Figure  A. 459.  Pd  vs  Time 


Figure  A. 460.  Pd  vs  Time 
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Figure  A. 464.  Pd  vs  Time 
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Figure  A. 465.  Pd  vs  Time 


Figure  A. 467.  Pd  vs  Time 
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Figure  A. 468.  Pd  vs  Time 
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Figure  A. 470.  Pd  vs  Time 
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Figure  A. 473.  Pd  vs  Time 
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Figure  A. 474.  Pd  vs  Time 
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Figure  A. 477.  Pd  vs  Time  Figure  A. 478.  Pd  vs  Time 


Figure  A. 479.  Pd  vs  Time 


Figure  A. 480.  Pd  vs  Time 
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APPENDIX  B: 
Tables 


B.l  Discrete 

Table  B.l.  Parameter  Variation  for  p,  for  Uniform  Distribution 


Parameter 

b 

r2 

SSE 

RMSE 

P  uni  form  ^uniform 

0.01061 

0.99242 

0.006608 

0.027096 

Puni  form  iQ  long 

0.009988 

0.99935 

0.00049844 

0.007442 

Puniform  ■>  corner 

0.01198 

0.99795 

0.0015533 

0.013137 

Best  fit  line  parameters  for  a  uniformly  distributed  location  probability  density  of 
the  form  f(x)  =  1  -  e~bx 


Table  B.2.  Parameter  Variation  for  pi  for  Corner  Stratification 


Parameter 

b 

r2 

SSE 

RMSE 

Pcornerstrat  ^uniform 

0.01321 

0.99908 

0.00065853 

0.008554 

Pcornerstrat  iQlong 

0.00896 

0.9992 

0.00066179 

0.0085751 

Pcornerstrat  ^corner 

0.01119 

0.99808 

0.0013515 

0.012254 

Best  fit  line  parameters  for  a  corner  stratified  distribution  location  probability 
density  of  the  form  f(x)  =  1  -  e~bx 
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Table  B.3.  Parameter  Variation  for  p,  for  Centered  Gaussian 


Parameter 

b 

r2 

SSE 

RMSE 

P  Gauss  ^uniform 

0.0147 

0.99781 

0.0014554 

0.012717 

PGaussiQlong 

0.01315 

0.99609 

0.0026167 

0.017051 

P Gauss  ^corner 

0.01589 

0.9963 

0.00234 

0.016124 

Best  fit  line  parameters  for  a  centered  Gaussian  distribution  location  probability 
density  of  the  form  f(x)  =  1  -  e~bx 


Table  B.4.  Parameter  Variation  for  pi  for  Side  Gaussian 


Parameter 

b 

r2 

SSE 

RMSE 

P  Gauss  ^uniform 

0.02405 

0.9944 

0.0027756 

0.017561 

P  GaussiQlong 

0.01428 

0.99225 

0.0047039 

0.0228621 

P Gauss  ^corner 

0.01895 

0.98981 

0.0054218 

0.024544 

Best  fit  line  parameters  for  a  offset  Gaussian  distribution  location  probability  den¬ 
sity  of  the  form  f(x)  =  1  -  e~bx 


B.2  Bayesian 

Table  B.5.  Parameter  Variation  for  p ,•  for  Uniform  Distribution 


Parameter 

a 

b 

r2 

SSE 

RMSE 

P  strut  ^uniform 

0.90643 

-0.01139 

0.99645 

0.004106 

0.015103 

P  strati  long 

0.89218 

-0.011534 

0.98964 

0.0104 

0.024037 

P  st  rat  corner 

0.8876 

-0.014151 

0.99384 

0.0058881 

0.018086 

Best  fit  line  parameters  for  a  uniformly  distributed  location  probability  density  of 
the  form  f(x)  —  a  -  e~bx 
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Table  B.6.  Parameter  Variation  for  p,  for  Corner  Stratification 


Parameter 

a 

b 

r2 

SSE 

RMSE 

Pstratitfwiiform 

0.88918 

-0.018268 

0.976 

0.0175 

0.031181 

Pstrat  ?Qlong 

0.89054 

-0.012323 

0.98204 

0.016665 

0.030427 

Pstrat  ^corner 

0.89926 

-0.011051 

0.99147 

0.008804 

0.022116 

Best  fit  line  parameters  for  a  corner  stratified  distribution  location  probability 
density  of  the  form  f(x)  =  a  -  e~bx 


Table  B.7.  Parameter  Variation  for  pi  for  Centered  Gaussian 


Parameter 

a 

b 

r2 

SSE 

RMSE 

P  strut  5  ^  uniform 

0.88971 

-0.017384 

0.9847 

0.012037 

0.025859 

PstrahQlong 

0.88274 

-0.016095 

0.98036 

0.01587 

0.029693 

Pstrat  corner 

0.88519 

-0.019317 

0.98268 

0.01265 

0.02651 

Best  fit  line  parameters  for  a  centered  Gaussian  distribution  location  probability 
density  of  the  form  f(x)  =  a  -  e~bx 


Table  B.8.  Parameter  Variation  for  pi  for  Side  Gaussian 


Parameter 

a 

b 

r2 

SSE 

RMSE 

Pstrat  ^uniform 

0.88957 

-0.029116 

0.97623 

0.012426 

0.026275 

PstratiQ  long 

0.88243 

-0.01768 

0.96822 

0.023033 

0.035772 

Pstrat  corner 

0.88301 

-0.023588 

0.96697 

0.019693 

0.0330771 

Best  fit  line  parameters  for  a  offset  Gaussian  distribution  location  probability  den¬ 
sity  of  the  form  f(x)  =  a  -  e~bx 
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B.3  Continuous 

B.3.1  Inline  Ladder  Search 

Table  B.9.  Inline  Ladder  Search:  Time  to  Po  =  90%(V  =  5) 


Velocities 

8 

9 

Ship  1 

10 

11 

12 

8 

14.8081 

14.5798 

15.6685 

16.1428 

16.6395 

9 

15.7111 

15.6258 

16.6558 

16.8813 

17.6232 

Ship  2  10 

18.5017 

18.0867 

17.8148 

18.2146 

18.7484 

11 

18.2434 

19.5095 

19.5958 

18.8285 

20.0086 

12 

18.1569 

19.1277 

20.5486 

21.9053 

21.3588 

Inline  Ladder  Search 


Table  B.10.  Ini 

ine  Ladder  Search: 

Time  to  Pd 

=  90%  (V 

=  6) 

Ship  1 

Velocities 

8 

9 

10 

11 

12 

8 

12.7567 

12.2523 

12.1515 

12.4612 

13.1644 

9 

13.3414 

12.6024 

12.3187 

12.5611 

13.1930 

Ship  2  10 

14.0466 

13.0154 

12.6810 

13.2800 

13.5936 

11 

13.7793 

15.2770 

14.1587 

14.0189 

14.6239 

12 

13.9657 

15.0961 

16.4268 

15.6398 

15.5540 

Inline  Ladder  Search 
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Table  B.ll.  Inline  Ladder  Search:  Time  to  Pd  =  90%(V  =  7) 


Velocities 

8 

9 

Ship  1 

10 

11 

12 

8 

10.9647 

10.6280 

10.4097 

10.2503 

10.2367 

9 

11.7028 

11.2848 

10.7237 

10.5538 

10.5661 

Ship  2  10 

12.7877 

11.5763 

11.0378 

10.8235 

10.7182 

11 

12.2397 

12.4420 

11.3968 

11.0259 

10.9716 

12 

11.8331 

12.0358 

12.8213 

11.4023 

11.4783 

Inline  Ladder  Search 


Table  B.12.  Inline  Ladder  Search:  Time  to  Pd  =  90%(V  =  8) 


Velocities 

8 

9 

Ship  1 

10 

11 

12 

8 

8.7756 

8.6460 

8.8498 

8.7901 

8.8294 

9 

10.1282 

9.7988 

9.2987 

9.2396 

9.0239 

Ship  2  10 

11.4149 

10.2606 

9.8759 

9.5892 

9.3740 

11 

11.1874 

11.3872 

10.2884 

9.8888 

9.6073 

12 

10.9687 

11.0076 

11.0854 

10.1519 

9.7240 

Inline  Ladder  Search 
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Table  B.13.  Inline  Ladder  Search:  Time  to  Pd  =  90%(V  =  9) 


Velocities 

8 

9 

Ship  1 

10 

11 

12 

8 

6.8662 

6.9782 

7.1117 

7.5250 

7.7083 

9 

7.5069 

7.8241 

7.7043 

7.8639 

7.9255 

Ship  2  10 

9.7362 

9.0589 

8.7202 

8.4741 

8.3054 

11 

10.1011 

10.0333 

9.2198 

8.9486 

8.6854 

12 

10.0425 

10.1511 

10.0968 

9.2826 

8.9569 

Inline  Ladder  Search 


Table  B.14.  Inline  Ladder  Search:  Time  to  Pd  =  99%(V  =  10) 


Velocities 

8 

9 

Ship  1 

10 

11 

12 

8 

6.6247 

6.4570 

6.3438 

6.4058 

6.5204 

9 

6.5861 

6.3278 

6.3438 

6.4924 

6.8038 

Ship  2  10 

7.2954 

6.9782 

7.0840 

7.0118 

7.1582 

11 

9.1759 

8.6565 

8.2237 

7.8774 

7.6543 

12 

9.2135 

9.2844 

9.0009 

8.3631 

8.1504 

Inline  Ladder  Search 


176 


B.3.2  Multipath  Ladder 


B.3.3  Multipath  Ladder  Search 

Table  B.15.  Multipath  Ladder  Search:  Time  to  Pd  =  90% (V  =  5) 


Velocities 

8 

9 

Ship  1 

10 

11 

12 

8 

14.6171 

14.9854 

15.6398 

16.1907 

16.5752 

9 

14.9427 

15.8392 

16.3410 

16.9964 

17.6940 

Ship  2  10 

18.2155 

17.4714 

17.5716 

18.2338 

18.7098 

11 

18.2913 

19.0874 

18.5983 

18.9148 

19.9443 

12 

18.3240 

19.1792 

20.6129 

20.6580 

21.1080 

Multipath  Ladder  Search 


Table  B.16.  Multipath  Ladder  Search:  Time  to  Pd  =  90%(V  =  6) 


Velocities 

8 

9 

Ship  1 

10 

11 

12 

8 

12.7024 

12.2912 

12.2072 

13.0803 

13.1644 

9 

12.9135 

12.4857 

12.4580 

12.9206 

13.2073 

Ship  2  10 

13.7122 

12.7089 

12.7089 

13.2600 

13.7510 

11 

14.0389 

14.0988 

13.6395 

14.2585 

14.5953 

12 

13.9371 

15.6255 

15.8545 

15.2106 

15.5254 

Multipath  Ladder  Search 
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Table  B.17.  Multipath  Ladder  Search:  Time  to  Pd  =  90°/c(V  =  7) 


Velocities 

8 

9 

Ship  1 

10 

11 

12 

8 

11.4414 

10.9266 

10.4994 

10.2503 

10.3888 

9 

11.8222 

11.1654 

10.7686 

10.6212 

10.7942 

Ship  2  10 

12.3390 

11.3070 

10.9481 

10.8235 

11.0729 

11 

12.4083 

11.7002 

11.1270 

10.8910 

11.3263 

12 

12.0611 

12.3145 

11.5037 

11.0729 

11.630 

Multipath  Ladder  Search 


Table  B.18.  Multipath  Ladder  Search:  Time  to  Pd  =  90%(V  =  8) 


Velocities 

8 

9 

Ship  1 

10 

11 

12 

8 

10.1501 

9.7165 

9.2987 

9.0398 

8.9850 

9 

10.7870 

10.2105 

9.7476 

9.2895 

9.1406 

Ship  2  10 

11.2226 

10.5813 

9.9400 

9.5392 

9.45180 

11 

11.5370 

10.7379 

10.0886 

9.7890 

9.6073 

12 

11.0854 

11.2021 

10.3464 

9.8407 

9.6462 

Multipath  Ladder  Search 
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Table  B.19.  Multipath  Ladder  Search:  Time  to  Pd  =  90%(V  =  9) 


Velocities 

8 

9 

Ship  1 

10 

11 

12 

8 

8.1866 

8.2470 

8.2123 

8.1351 

8.0883 

9 

9.3043 

8.9871 

8.7202 

8.3385 

8.1969 

Ship  2  10 

10.3288 

9.5669 

9.2282 

8.7452 

8.4140 

11 

10.7790 

10.0333 

9.4231 

9.0164 

8.6854 

12 

10.4768 

10.5311 

9.7168 

9.1197 

8.7940 

Multipath  Ladder  Search 


Table  B.20.  Multipath  Ladder  Search:  Time  to  Pd  =  90%(V  =  10) 


Velocities 

8 

9 

Ship  1 

10 

11 

12 

8 

6.7825 

6.8444 

6.7668 

7.0983 

7.0873 

9 

6.9735 

7.4901 

7.4011 

7.5312 

7.3708 

Ship  2  10 

8.7756 

8.4585 

8.2470 

7.8774 

7.6543 

11 

9.9550 

9.1759 

8.6565 

8.3102 

8.0087 

12 

9.8514 

9.8514 

9.0009 

8.5757 

8.221 

Multipath  Ladder  Search 
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B.3.4 


Spiral 


B.3.5  Spiral  Search 

Table  B.21.  Spiral  Search:  Time  to  Po  =  90%(V  =  5) 


Velocities 

8 

9 

Ship  1 

10 

11 

12 

8 

19.0117 

19.4895 

19.5176 

19.7109 

20.2015 

9 

21.6669 

20.4074 

20.8340 

20.8907 

20.0729 

Ship  2  10 

22.2220 

23.1092 

21.5638 

21.5909 

21.67388 

11 

21.8115 

22.8090 

25.0000 

22.5980 

22.65114 

12 

21.7960 

22.4003 

23.4676 

25.0000 

25.0000 

Spiral  Search 


Table  B.22.  Spiral  Search:  Time  to  Pd  =  90%  (V  =  6) 


Velocities 

8 

9 

Ship  1 

10 

11 

12 

8 

13.8424 

14.3527 

14.3253 

14.5182 

16.0405 

9 

16.0252 

15.7918 

16.2205 

16.5751 

15.0818 

Ship  2  8 

18.3108 

17.6140 

17.1124 

17.3739 

17.5859 

11 

18.2326 

18.7518 

18.6519 

18.2126 

18.1725 

12 

18.3013 

18.9309 

19.5462 

19.8467 

18.9309 

Spiral  Search 
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Table  B.23.  Spiral  Search:  Time  to  Pd  =  90%(V  =  7) 


Velocities 

8 

9 

Ship  1 

10 

11 

12 

8 

11.5209 

11.8819 

12.1147 

12.1723 

12.5932 

9 

11.7625 

12.0013 

12.3839 

12.1723 

11.4023 

Ship  2  10 

15.1658 

13.7749 

13.5505 

13.7233 

14.0122 

11 

15.0720 

15.8476 

15.0383 

14.8697 

15.1017 

12 

14.8737 

16.1153 

16.5460 

15.9886 

15.7099 

Spiral  Search 


Table  B.24.  Spiral  Search:  Time  to  Pd  -  90%(V  =  8) 


Velocities 

8 

9 

Ship  1 

10 

11 

12 

8 

10.3616 

10.2929 

10.9019 

11.0875 

11.2021 

9 

10.4576 

10.2929 

10.6454 

10.6879 

10.5409 

Ship  2  10 

11.8638 

10.5813 

10.6454 

10.7878 

10.6575 

11 

12.2362 

13.2850 

11.9365 

11.8866 

11.7855 

12 

12.0189 

13.6137 

13.9637 

13.1469 

12.9524 

Spiral  Search 
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Table  B.25.  Spiral  Search:  Time  to  Pd  =  90%(V  =  9) 


Velocities 

8 

9 

Ship  1 

10 

11 

12 

8 

9.6390 

9.3043 

9.6515 

10.2367 

10.3139 

9 

9.8330 

9.4100 

9.3129 

9.7621 

9.8254 

Ship  2  10 

10.4135 

9.4822 

9.3129 

9.5587 

9.6625 

11 

10.5078 

10.0333 

9.5587 

9.6943 

9.71684 

12 

10.3682 

11.6710 

11.5625 

10.6397 

10.5854 

Spiral  Search 


Table  B.26.  Spiral  Search:  Time  to  Pd  =  90%(V  =  10) 


Velocities 

8 

9 

Ship  1 

10 

11 

12 

8 

9.3062 

8.6524 

8.6699 

9.5221 

9.4262 

9 

9.2981 

8.7815 

8.5642 

8.8296 

9.2135 

Ship  2  10 

9.6215 

8.8814 

8.5642 

8.5699 

8.8592 

11 

9.6087 

9.1759 

8.5699 

8.6565 

8.7883 

12 

9.3553 

9.7805 

9.0009 

8.7174 

8.788 

Spiral  Search 
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APPENDIX  C: 
Codes 


C.l 

Discrete 

C.1.1 

Main  Code 

i  clc 

2  clear 

3  discrete_search_prob 

4  discrete_loc_probs2 

5  close  all 

6  %p-uniform 

7  dis_pUni_qUni 

s  dis_pUni_qlongstrat 
s  dis_pUni_qbothstrat 

10  %p-long  stratified 

11  dis_plong  strat_qUni 

12  dis_plongstrat_qlongstrat 

13  dis_plongstrat_qbothstrat 

14  %p-both  stratified 

is  dis_pbothstrat_qUni 
i6  dis_pbothstrat_qlongstrat 
n  dis_pbothstrat_qbothstrat 
is  %p-normal  center 

19  dis_pnormwide_qUni 

20  dis_normwide_qlongstrat 

21  dis_normwide_qbothstrat 

22  %p- normal  side 

23  dis_pnormside_qUni 

24  dis_normside_qlongstrat 

25  dis_normside_qbothstrat 
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26  %p-normal  corner 

27  dis_pnormc_qUni 

28  dis_normc_qlongstrat 

29  dis_normc_qbothstrat 

30  m> 

31 

32  pd  =  [  Pd_uni_uni  ( end  )  ,  Pd_uni_ls  (  end  )  , Pd_uni_bs ( end )  ,Pd_ls_uni( 

end )  , . .  . 

33  Pd_ls_ls  ( end  )  ,  Pd_ls_b s  (  end  )  ,  Pd_bs_uni  ( end  )  ,  Pd_b s_ls  (  end  ) 

, Pd_bs_bs ( end )  . .  . 

34  ,  Pd_nw_uni  ( end  )  ,  Pd_nw_ls  (  end  )  ,  Pd_nw_bs  (  end  )  ,  Pd_ns_uni  ( 

end  )  , .  .  . 

35  Pd_ns_ls  ( end  )  ,  Pd_ns_bs  (  end  )  ,  Pd_nc_uni  ( end  )  ,  Pd_nc_ls  (  end  ) 

5  •  •  * 

36  Pd_nc_bs  ( end  )  ] 

37 

38  m 

39  %  close  all 

40  set(0,’defaulttextinterpreter’  , ’latex') 

41  clear  g 

42  [  ~  ,  Pd_uni_uni  ,  k]  =  discrete_search2  (  P_uni  ,  q_uni  ,500  ,  1) 

43  [  ~  ,  Pd_uni_ls  ,k]  =  discrete_search2  (  P_uni  ,  q_long_strat 

,500,  1  ); 

44  [  ~  ,  Pd_uni_bs  ,k]  =  discrete_search2  (  P_uni  ,  q_both_strat 

,500,  1  ); 

45  figure 

46 

47  stepl  =10:  round(. 05*  length  (Pd_uni_uni))  :  length(Pd_uni_uni)  ; 

48  step2  =  10:  round  (.05*  length  (Pd_uni_ls))  :  length(Pd_uni_ls)  ; 

49  step3  =  10:round(.05*length(Pd_uni_bs))  :  length(Pd_uni_bs)  ; 

50  plot(stepl  ,Pd_uni_uni(stepl  )  ,  ’  xk  ’  , step  2  ,Pd_uni_ls(step2)  ,  ’  xr 

,  step  3  ,  Pd_uni_bs  (step3)  ,  ’  xg  ’  ) 
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51  g  =  fittype  (’  1-exp  (b*x)  dependent  Prob  independent  ’x 

52 

53  [  fit_bb  1  ,  gof_bb  1]=  fit  (stepl  Pd_uni_uni  (stepl  )  ’  ,g,  ’ 

StartPoint  ’  ,[0]  ,  ’  Upper  ’  ,  [0] )  ; 

54  [  fit_bb2  ,  gof_bb2]=  fit(step2  ’  ,Pd_uni_ls(step2)  ’  ,g,  ’StartPoint 

’  ,  [0]  ,  ’Upper’  ,[0])  ; 

55  [  f  i  t_b  b  3  ,  gof_bb3  ]  =  fit  (step3  ’  ,  Pd_uni_bs  (step3  )  ’  ,g,  ’  StartPoint 

’  ,[0]  ,  ’Upper’  ,  [0] )  ; 

56  hold  on 

57  plot  (  fit_bb  1  ,  ’ k  ’ ) 

58  plot  (  fit_bb2  ,  ’  r  ’ ) 

59  plot  (  f i t_b b 3  ,  ’  g  ’ ) 
so  xlabel(’n’) 

61  ylabel  (  ’  $P_{D}$  ’ ) 

62  t  i  1 1  e  (  ’  $P_ {D }  $  vs  Number  of  Glimpses  and  Exponential  Fitted 

Curves  ’ ) 

63  legend  (’ (p_{  uniform  }  ,q_{  uniform  })’,’ (p_{  uniform  }  ,q_{  s  tr  at  a 

})  ’  ,  ’  (p_{  uniform  }  ,q_{  str at  b  })  ’  ... 

64  ,  sprintf('l-eA{%.5f*x}’  ,  f  i  t  _  b  b  1  .b)  ,  sprintf(’l-eA{c7o.5f*x} 

’  ,  fit_bb2  .b)  ,... 

65  sprintf(’l-eA{%.5f*x}’,fit_bb3.b),  ’location’,  ’best’) 

66  print  (  ’  Bayes_Uni_Fit  ’  ,  ’-dpng  ’ ) 

67 

68  [  ~  ,  Pd_ls_uni  ,  k]  =  discrete_search2  (  P_long_strat  ,  q_uni 

,500,  1  ); 

69  [  ~  ,  Pd_ls_ls  ,k]  =  discrete_search2  (  P_long_strat  , 

q_long_strat  ,500 ,  1  ); 

70  [  ~  ,~  ,~  ,  Pd_ls_bs  ,k]  =  discrete_search2  (  P_long_strat  , 

q_both_strat  ,500 ,  1  ); 

71  figure 

72  stepl=10:round(.05*length(Pd_ls_uni))  :  length(Pd_ls_uni)  ; 

73  step2  =  10:  round  (.05*  length  (Pd_ls_ls))  :  length(Pd_ls_ls)  ; 
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74  step3  =  10:round(.05*length(  Pd_ls_bs  )):  length  (  Pd_ls_bs  )  ; 

75  plot(stepl  ,  Pd_ls_uni(stepl)  ,’xk’  , step 2  ,Pd_ls_ls(step2)  ,  ’  xr  ’  , 

step 3  ,  Pd_ls_bs  (  step 3  )  ,  ’  xg  ’ ) 

76  [  fi  t_bb4  ,  gof_bb4  ]=  fit  (stepl  ’  ,Pd_ls_uni(stepl  )  ’  ,g,  ’  StartPoint 

’  ,[0]  ,  ’Upper’  ,[0])  ; 

77  [  f  i  t_b  b  5  ,  gof_bb5  ]=  fit  (step2  ’  ,Pd_ls_ls(step2)  ’  ,g  ,  ’StartPoint  ’ 

,[0]  ,  ’Upper’  ,[0])  ; 

78  [  fit_bb6  ,  gof_bb6  ]=fit(step3  ’  ,Pd_ls_bs(step3)  ’  ,g  ,  ’StartPoint  ’ 

,[0]  ,  ’Upper’  ,[0])  ; 

79  hold  on 

so  plot  (  fit_bb4  ,  ’ k  ’ ) 
si  plot  (  f i t_b b 5  ,  ’  r  ’ ) 

82  plot  (  fit_bb6  ,  ’  g  ’ ) 

83  xlabel(’n’) 

84  ylabel  (  ’  $P_{D}$  ’ ) 

85  t  i  1 1  e  (  ’  $P_ {D }  $  vs  Number  of  Glimpses  and  Exponential  Fitted 

Curves  ’ ) 

86  legend  (’ (p_{  s  tr  at  a  },  q_{  uniform  })’,’( p_  {  s  t  r  a  t  a},q_{strat  a 

}) ’  , ’ (p_{ strat  a } , q_ { strat  b  })  ’  .  .  . 

87  ,  sprintf('l-eA{%.5f*x}’  ,fit_bb4.b)  ,  sprintf(’l-eA{c7o.5f*x} 

’  ,fit_bb5  .b) 

88  sprintf('l-eA{%.5f*x}’  ,fit_bb6.b)  ,  ’location’  ,  ’  best  ’  ) 
so  print  (  ’  Bayes_Long_Fit  ’  ,  ’-dpng  ’ ) 

90 

91  [  ~  ,  Pd_bs_uni  ,  k]  =  discrete_search2  (  P_both_strat  ,  q_uni 

,500,  1  ); 

92  [  ~  ,  Pd_bs_ls  ,k]  =  discrete_search2  (  P_both_strat  , 

q_long_strat  ,500 ,  1  ); 

93  [  ~  ,~  ,~  ,  Pd_bs_bs  ,k]  =  discrete_search2  (  P_both_strat  , 

q_both_strat  ,500 ,  1  ); 

94  figure 

95  step  1  =10:  round(. 05*  length  (Pd_bs_uni))  :  length(Pd_bs_uni)  ; 

96  step2  =  10:round(.05*length(Pd_bs_ls))  :  length(Pd_bs_ls)  ; 
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97  step3  =  10:round(.05*length(  Pd_bs_bs  )):  length  (  Pd_bs_bs  )  ; 

98  plot(stepl,  Pd_bs_uni  (stepl)  ,’xk’  , step  2  ,  Pd_bs_ls  (step2)  ,  ’  xr  ’  , 

step  3  ,  Pd_bs_bs  (step3)  ,’xg’) 

99  [  fit_bb7  ,  gof_bb7  ]=  fit  (stepl  Pd_bs_uni  (stepl  )  ’  ,g,  ’  StartPoint 

’  ,[0]  ,  ’Upper’  ,[0])  ; 

100  [  fit_bb8  ,  gof_bb8  ]=fit(step2  ’  ,Pd_bs_ls(step2)  ’  ,g  ,  ’StartPoint  ’ 

,[0]  ,  ’Upper’  ,[0])  ; 

101  [  fit_bb9  ,  gof_bb9  ]=fit(step3’,  Pd_bs_bs  (step3)  ’  ,g,  ’  StartPoint  ’ 

,[0]  ,  ’Upper’  ,[0])  ; 

102  hold  on 

103  plot(fit_bb7  ,  ’k’) 

104  plot  (  fit_bb8  ,  ’  r  ’ ) 
los  plot  (  fit_bb9  ,  ’  g  ’  ) 
ioo  x label  (  ’n  ’  ) 

107  ylabel  (  ’  $P_{D}$  ’ ) 

ms  title  (  ’  $P_ { D }  $  vs  Number  of  Glimpses  and  Exponential  Fitted 
Curves  ’ ) 

109  legend  (’ (p_{  s  tr  at  b  },  q_{  uniform  })’,’( p_  {  s  t  r  a  t  b},q_{strat  a 
})  ’  ,  ’  (p_{  strat  b  }  ,  q_  {  s  t  r  a  t  b})’... 
no  ,  sprintf(’l-eA{%.5f*x}’  ,  fit_bb  7  .b)  ,  sprintf('l-eA{%.5f*x} 

’  ,  fit_bb  8  .b) 

in  sprintf(’l-eA{%.5f*x}’  ,fit_bb9.b)  ,  ’location’  ,  ’  best  ’  ) 

n2  print  (  ’Bayes_Both_Fit  ’  ,  ’-dpng  ’ ) 

113 

n4  [  ~  , ~  , Pd_nw_uni  ,  k]  =  discrete_search2  (P_norm_wide  ,  q_uni 

,500,  1  ); 

ns  [  ~  ,  Pd_nw_ls  ,k]  =  discrete_search2  (P_norm_wide  , 

q_long_strat  ,500 ,  1  ); 

lie  [  ~  ,~  ,~  , Pd_nw_bs  ,k]  =  discrete_search2  (P_norm_wide  , 
q_both_strat  ,500  ,  1  ); 
in  figure 

ns  stepl=10:round(.05*length(  Pd_nw_uni  )):length(  Pd_nw_uni )  ; 
no  step2  =  10:round(.05*length(  Pd_nw_ls  )):length(  Pd_nw_ls  )  ; 
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120  step3  =  10:  round(. 05*  length  (Pd_nw_bs) )  :  length  (Pd_nw_bs)  ; 
m  plot  (  step  1  ,  Pd_nw_uni  (stepl)  ,’xk’  , step  2  ,  Pd_nw_ls  (step2)  ,  ’  xr  '  , 
step  3  ,  Pd_nw_bs  (step3)  ,’xg’) 

122  [  fit_bb  1  0  ,  gof_bb 10  ]  =  f  i  t  (  step  1  ’  ,  Pd_nw_uni  (stepl)  ’  ,g  , 

StartPoint  ’  ,[0]  ,  ’Upper’  ,  [  0  ] )  ; 

123  [  fit_bb  1  1  ,  gof_bbl  1  ]  =  f  it  (  step 2  ’  ,  Pd_nw_ls  (  step 2  )  ’  ,g  , 

StartPoint  ’  ,[0]  ,  ’  Upper  ’  ,  [0] )  ; 

124  [  fit_bb  1  2  ,  gof_bb  12]  =  fit  (step  3  ’  ,  Pd_nw_bs  (step3)  ’  ,g  , 

StartPoint  ’  ,[0]  ,  ’  Upper  ’  ,  [0] )  ; 

125  hold  on 

126  plot  (  fit_bb  1  0  ,  ’k’) 

127  plot  (  fit_bb  1  1  ,’r’) 

128  plot  (  fit_bb  1  2  ,  ’g’) 
m  x label ( ’n ’ ) 

no  ylabel  (  ’  $P_{D}$  ’ ) 

131  title  (  ’  $P_ { D }  $  vs  Number  of  Glimpses  and  Exponential  Fitted 

Curves  ’ ) 

132  legend  (’  (p_{norm  }  ,q_{  uniform  })’,’  (p_{norm  }  ,q_{  strat  a})’,’( 

p_  {  norm  }  ,  q_  {  strat  b  } )  ’  .  .  . 

133  ,sprintf(’l-eA{%.5f*x}’  ,fit_bbl0.b)  ,  sprintf(’l-eA{%.5f*x 

}  ’  ,  fit_bb  1  1  .b)  , .  .  . 

134  sprintf(’l-eA{%.5f*x}’  ,fit_bbl2.b)  ,  ’location’  ,  ’  bes  t  ’  ) 

135  print  (  ’  Bayes_Norm_Fit  ’  ,  ’-dpng  ’ ) 

136 

137  [  ~  ,  Pd_ns_uni  ,  k]  =  discrete_search2  (  P_norm_side  ,  q_uni 

,500,  1  ); 

ns  [  ~  ,  Pd_ns_ls  ,k]  =  discrete_search2  (  P_norm_side  , 

q_long_strat  ,500 ,  1  ); 

139  [  ~  ~ , Pd_ns_bs  ,k]  =  discrete_search2  (  P_norm_side  , 

q_both_strat  ,500 ,  1  ); 
wo  figure 

mi  stepl=10:round(.05*length(Pd_ns_uni))  :  length(Pd_ns_uni)  ; 

W2  step2  =  10:round(.05*length(Pd_ns_ls))  :  length(Pd_ns_ls)  ; 
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143  step3  =  10:round(.05*length(  Pd_ns_bs  )):length(  Pd_ns_bs  )  ; 

144  plot  (  step  1  ,  Pd_ns_uni  (stepl)  ,’xk’  , step  2  ,Pd_ns_ls(step2)  ,  ’  xr  '  , 

step  3  ,  Pd_ns_bs  (step3)  ,’xg’) 

ms  [  fit_bb  1  3  ,  gof_bb  13]  =  fit  (step  1  Pd_ns_uni  (stepl)  ’  ,g  , 

StartPoint  ’  ,[0]  ,  ’Upper’  ,[0])  ; 
we  [  fit_bb  1  4  ,  gof_bb  14  ]=  fit  (  step 2  ’  ,  Pd_ns_ls  (  step 2  )  ’  ,g  , 

StartPoint  ’  ,[0]  ,  ’  Upper  ’  ,  [0] )  ; 

147  [  fit_bb  1  5  ,  gof_bb  15]  =  fit  (step  3  ’  ,  Pd_ns_bs  (step3)  ’  ,g  , 

StartPoint  ’  ,[0]  ,  ’  Upper  ’  ,  [0] )  ; 

148  hold  on 

149  plot  (  fit_bb  1  3  ,  ’k’) 
iso  plot(fit_bbl4,’r’) 

151  plot  (  fit_bb  15  ,  ’ g  ’ ) 

152  x label  (  ’n  ’  ) 

153  ylabel  (  ’  $P_{D}$  ’ ) 

154  title  (  ’  $P_ { D }  $  vs  Number  of  Glimpses  and  Exponential  Fitted 

Curves  ’ ) 

155  legend  (’ (p_{  norm  b  },  q_{  uniform  })’,’( p_  {  norm  b],q_{strat  a})’ 

,’(p_{norm  b},q_{strat  b})’... 

156  ,sprintf(’l-eA{%.5f*x}’  ,  fit_bb  1  3  .b)  ,  sprintf(’l-eA{%.5f*x 

}  ’  ,  fit_bb  1 4  .b)  , . . . 

157  sprintf(’l-eA{%.5f*x}’  ,  fit_bb  1  5  .b)  ,  ’location’  ,  ’  bes  t  ’  ) 

158  print  (  ’  Bayes_NormSide_Fit  ’  ,  ’-dpng  ’ ) 

159 

wo  [  ~  ,  Pd_nc_uni  ,  k]  =  discrete_search2  (  P_norm_corner  ,  q_uni 

,500,  1  ); 

161  [  ~  ,  Pd_nc_ls  ,k]  =  discrete_search2  (  P_norm_corner  , 

q_long_strat  ,500 ,  1  ); 

162  [  ~  ~  ,  Pd_nc_bs  ,k]  =  discrete_search2  (  P_norm_corner  , 

q_both_strat  ,500 ,  1  ); 

163  figure 

164  stepl=10:round(.05*length(Pd_nc_uni))  :  length(Pd_nc_uni)  ; 

165  step2  =  10:round(.05*length(Pd_nc_ls))  :  length(Pd_nc_ls)  ; 
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166  step3  =  10:round(.05*length(  Pd_nc_bs  )):length(  Pd_nc_bs  ) 

167  plot  (  step  1  ,  Pd_nc_uni  (stepl)  ,’xk’  , step  2  ,Pd_nc_ls(step2)  ,  ’  xr  '  , 

step  3  ,  Pd_nc_bs  (step3)  ,’xg’) 

168  [  fit_bb  1  6  ,  gof_bb  16]  =  fit  (step  1  Pd_nc_uni  (stepl)  ’  ,g  , 

StartPoint  ’  ,[0]  ,  ’Upper’  ,  [  0  ] )  ; 

169  [  fit_bb  1  7  ,  gof_bb  17]  =  fit  (step  2  ’  ,Pd_nc_ls(step2)  ’  ,g  , 

StartPoint  ’  ,[0]  ,  ’  Upper  ’  ,  [0] )  ; 
no  [  fit_bb  1  8  ,  gof_bb  18]  =  fit  (step  3  ’  ,  Pd_nc_bs  (step3)  ’  ,g  , 

StartPoint  ’  ,[0]  ,  ’  Upper  ’  ,  [0] )  ; 

171  hold  on 

172  plot  (  fit_bb  1  6  ,  ’k’) 

173  plot(fit_bbl7,’r’) 
m  plot  (  fit_bb  18  ,  ’ g  ’ ) 
ns  x label  (  ’n  ’  ) 

n6  ylabel  (  ’  $P_{D}$  ’ ) 

177  title  (  ’  $P_ { D }  $  vs  Number  of  Glimpses  and  Exponential  Fitted 
Curves  ’ ) 

ns  legend  (’ (p_{  norm  c  },  q_{  uniform  })’,’( p_  {  norm  c],q_{strat  a})’ 
,’(p_{norm  c},q_{strat  b})’... 

no  ,sprintf(’l-eA{%.5f*x}’  ,fit_bbl6.b)  ,  sprintf(’l-eA{%.5f*x 

}  ’  ,  fit_bb  1  7  .b)  , . . . 

iso  sprintf(’l-eA{%.5f*x]’,fit_bbl8.b),  ’location’,  ’best’) 

lsi  print  (  ’  Bayes_NormCorner_Fit  ,  ' -dpng  ’  ) 

182  format  shortG 

183  b  =  [fit_bbl  .b,  fit_bb2  .b,  f i t _ b b 3  .b,  fit_bb4  .b,  f i t _ b b 5  .b,  fit_bb6 

.  b  ,  fit_bb7  .  b  ,  .  .  . 

184  fit_bb 8  . b  ,  fit_bb9  .b,  fit_bb  10  . b  ,  fit_bb  1 1  .b,  fit_bb  1 2  .b, 

fit_bb  1 3  . b  ,  fit_bb  1 4  .b  , .  .  . 

185  fit_bb  1  5  .b,fit_bbl6  .b,fit_bbl7  .b,fit_bbl8  .b] 

186  r2  =  [gof_bbl  .  rsquare  ,  gof_bb2  .  rsquare  ,  gof_bb3  .  rsquare  ,  gof_bb4  . 

rsquare  ,  .  .  . 

is?  gof_bb5  .  rsquare  ,  gof_bb6  .  rsquare  ,  gof_bb7  .  rsquare  ,  gof_bb8  . 

rsquare  ,  .  .  . 
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gof_bb9  .  r  square  ,  gof_bb  lO.rsquare  ,  gof_bb  11  .rsquare  , 
gof_bbl2  .  rsquare  ,  .  .  . 

gof_bbl3  .  rsquare  ,  gof_bbl4  .  rsquare  ,  gof_bbl5  .  rsquare  , 
gof_bbl6  .  rsquare  ,  .  .  . 
gof_bbl7  .  rsquare  ,  gof_bbl8  .  rsquare  ] 

C.1.2  Search  Probabilities  cjj 

format  short 
format  compact 
Area  =  1000; 
x=-5:5;y=x; 

[X,Y]  =  meshgrid(x,y)  ; 

m> 

q0  = .  3  ; 

q_uni=q0*ones (10,10) ; 
figure 

subplot  (1,3,1) 

imagesc([  -  4.5  ,4.5]  ,[  -  4.5  ,4.5]  ,  q_uni  ,  [0  , 1  ] ) 

set(gca,  ’YDir’,  ’normal') 

cb_q=colorbar ; 

title  (cb_q  ,  ’ q_{ i } ’ ) 

hold  on 

plot  (X,Y,  ’k’)  ; 
plot  (Y,X,  ’k’)  ; 
t  i  1 1  e  (  ’  Uniform  ’ ) 
xlabel  (  ' Long  ’ ) 
y label  (  ’  Lat  ’ ) 

m 

% 

q0  = . 07  ; 

qx=l : 10; 

qy  =  ones  (10,1)  ; 
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q_long_strat=qy*qx*q0*.9; 
subplot (1 ,3 ,2) 

imagesc([  -  4.5  ,4.5]  ,[  -  4.5  ,4.5]  ,q_long_strat  ,[0  ,max( 
q_long_strat  (:))]) 
set(gca,  ’YDir’,  ’normal') 
cb_q=colorbar ; 
title  (cb_q  ,  ’q_{ i } ’ ) 
hold  on 

plot  (X,Y,  ’k’)  ; 
plot  (Y,X,  ’k’)  ; 
t i 1 1 e  (  ’  S tr at .  a’) 
xlabel  (  ’ Long  ’ ) 
y label  (  ’  Lat  ’  ) 

m 

q0  = .  0  1 ; 
qx  =  ( 1 : 1 0)  ; 

qy  =  ( 1 : 10) 

q_both_strat=qy*qx*q0*.9; 
subplot  (1  ,3  ,3) 

imagesc([  -  4.5  ,4.5]  ,[  -  4.5  ,4.5]  ,q_both_strat  ,[0  ,max( 
q_both_strat  (:))]) 
set(gca,  ’YDir’,  ’normal') 
cb_q=colorbar ; 
title  (cb_q  ,  ’ q_{ i } ’ ) 
hold  on 

plot  (X,Y,  ’k’)  ; 
plot  (Y,X,  ’k’)  ; 
t  i  1 1  e  (  ’  S  tr  at  .  b’) 
xlabel  (  ' Long  ’ ) 
y  label  (  ’Lat’) 

m 

print  (  ’Discrete_Miss_Probablities  ’  ,  ’-dpng  ’ ) 
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C.1.3  Location  Probabilities  pi 

%  clear 

close  all 
clc 

format  short 
format  compact 
Area  =  1000; 

Recreate  the  grid  positions 

p0  =  0; 

x  =  -  5 : 5  ;  y=x  ; 

[X,Y]  =  meshgrid(x,y)  ; 

m> 

P_uni  =  ones(10,10)  ; 

P_uni=(l-p0)*P_uni/  sum(  P_uni  ( : )  )  ; 
figure 

subplot  (2,3,1) 

imagesc([  -  4.5  ,4.5]  ,[  -  4.5  ,4.5]  ,  P_uni  ,  [0  , 1  ] ) 

set(gca,  ’YDir’,  ’normal’) 

cb_p=colorbar ; 

title (cb_p  ,  ’ p_ { i } ’ ) 

hold  on 

plot  (X,Y,  ’k’)  ; 
plot  (Y,X,  ’k’)  ; 
t  i  1 1  e  (  ’  Uniform  ’ ) 
xlabel ( ’ Long .  ’  ) 
y label  (  ’  Lat  ’ ) 

m 

Px  -  1 : 1 0 ; 

Py=ones  (10,1)  ; 

P_long_strat=Py*Px ; 

P_long_strat=(l-pO)*P_long_strat /sum(  P_long_strat  (:)  )  ; 
subplot (2 ,3 ,2) 

imagesc([  -  4.5  ,4.5]  ,[  -  4.5  ,4.5]  ,P_long_strat  ,[0  ,max( 


193 


33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 

51 

52 

53 

54 

55 

56 

57 

58 

59 

60 

61 

62 

63 


P_long_strat  (:))]) 
set(gca,  ’  YDir’,  ’normal') 
cb_p=colorbar ; 
title  (cb_p  ,  ’ p_{ i } ’ ) 
hold  on 

plot  (X,Y,  ’k’)  ; 
plot  (Y,X,  ’k’)  ; 
title  (  ’  Stratified  a  ’  ) 
xlabel ( ' Long .  ’ ) 
y label  (  ’  Lat  ’  ) 

m 

Px  =  1 : 1 0 ; 

Py  =  (l:10) 

P_both_strat=Py*Px ; 

P_both_strat=(l-pO)*P_both_strat /sum(  P_both_strat  (:)  )  ; 
%  figure 
subplot  (2  ,3  ,3) 

imagesc([  -  4.5  ,4.5]  ,[  -  4.5  ,4.5]  ,P_both_strat  ,[0  ,max( 
P_both_strat  (:))]) 
set(gca,  ’YDir’,  ’normal’) 
cb_p=colorbar ; 
title (cb_p  ,  ’ p_ { i } ’ ) 
hold  on 

plot  (X,Y,  ’k’)  ; 
plot  (Y,X,  ’k’)  ; 
title  (’Stratified  b’) 
xlabel ( ’ Long .  ’ ) 
y label  (  ’  Lat  ’  ) 

m 

Px  =  normpdf  ( [ -4 . 5  :  1 :  4 . 5  ]  ,0  , 2)  ; 

Py  =  normpdf ([-4. 5:1:4. 5], 0,5)  ’; 

P_norm_wide=Py*Px ; 

P_norm_wide=(l  -pO )  *P_norm_wide  /  sum(  P_norm_wide  ( : )  )  ; 
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subplot  (2,3  ,4) 

imagesc([  -  4.5  ,4.5]  ,[  -  4.5  ,4.5]  ,  P_norm_wide  ,  [  0  ,  max(  P_norm_wide 

(:))]) 

set(gca,  ’YDir’,  ’normal’) 
cb_p=colorbar ; 
title  (cb_p  ,  ’ p_ { i } ’ ) 
hold  on 

plot  (X,Y,  ’k’)  ; 
plot  (Y,X,  ’k’)  ; 
t  i  1 1  e  (  ’  Normal  a  ’ ) 
xlabel ( ’ Long .  ’  ) 
y label  (  ’  Lat  ’ ) 

m 

Px  =  normpdf  ( [ -4 . 5  :  1 :  4 . 5  ]  ,  3  , 2)  ; 

Py  =  normpdf ([-4. 5:1:4. 5], 3, 2) 

P_norm_corner=Py*Px ; 

P_norm_corner =(1  -pO )  *  P_norm_corner  /  sum(  P_norm_corner  ( : )  )  ; 
subplot  (2  ,3  ,5) 

imagesc([  -  4.5  ,4.5]  ,[  -  4.5  ,4.5]  ,  P_norm_corner  ,  [0  ,  max( 
P_norm_corner  (:))]) 
set(gca,  ’YDir’,  ’normal’) 
cb_p=colorbar ; 
title (cb_p  ,  ’ p_ { i } ’ ) 
hold  on 

plot  (X,Y,  ’k’)  ; 
plot  (Y,X,  ’k’)  ; 
title  (  ’  Normal  b  ’ ) 
xlabel ( ’ Long .  ’  ) 
y label  (  ’Lat  ’ ) 

m 

Px  =  normpdf  ( [ -4 . 5  :  1 :  4 . 5  ]  ,  3  , 2)  ; 

Py  =  normpdf ([-4. 5:1:4. 5], 0,2) 

P_norm_side=Py*Px ; 
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P_norm_side=(l  -pO )  *P_norm_side  /  sum(  P_norm_side  ( : )  )  ; 
subplot (2 ,3 ,6) 

imagesc([  -  4.5  ,4.5]  ,[  -  4.5  ,4.5]  ,  P_norm_side  ,  [  0  ,  max(  P_norm_side 

(:))]) 

set(gca,  ’YDir’,  ’normal') 
cb_p=colorbar ; 
title (cb_p  ,  ’ p_{ i } ’ ) 
hold  on 

plot  (X,Y,  ’k’)  ; 
plot  (Y,X,  ’k’)  ; 
t  i  1 1  e  (  ’  Normal  c  ’ ) 
xlabel ( ’ Long .  ’ ) 
y label  (  ’  Lat  ’  ) 

m 

print  (’Discrete_Location_Probablities  ’  ,’-dpng’) 

C.1.4  Sample  code  for  the  Uniform/Uniform  case 

[  pc  ,  n  ,  imax  ,  Pd_uni_uni  ,k]  =  di  s ere t e_s e arch 2  (  P_uni  ,  q_uni 
,50,  .95  ); 
figure 

subplot  (2,2,1) 

imagesc  ( [  -  4. 5, 4. 5], [-4. 5, 4. 5],  pc, [0  ,  max  (pc  (:))]) 

set(gca,  ’YDir’,  ’normal') 

cb_p=colorbar ; 

title (cb_p  ,  ’ p_ { i } ’ ) 

hold  on 

plot  (X,Y,  ’k’)  ; 
plot  (Y,X,  ’k’)  ; 
t  i  1 1  e  (  ’  Remainder  Prob  .  '  ) 
xlabel ( ' Long .  ’ ) 
y label  (  ’  Lat  .  ’ ) 

subplot (2 ,2 ,2) 
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imagesc([  -  4.5  ,4.5]  ,[  -  4.5  ,4.5]  ,n,[0  ,  max(n  (:))]) 

set(gca,  ’YDir’,  ’normal’) 

cb_p=colorbar ; 

title  ( cb_p  ,  ’N’ ) 

hold  on 

plot  (X,Y,  ’k’)  ; 
plot  (Y,X,  ’k’)  ; 
t  i  1 1  e  (’ Actual  Search’) 
xlabel ( " Long .  ’ ) 
y label  (  ’  Lat  .  ’  ) 
subplot (2 ,2 ,3) 

imagesc([  -  4.5  ,4.5]  ,[  -  4.5  ,4.5]  ,  P_uni  ,  [0  ,  max(  P_uni  ( : )  )  ] ) 

set(gca,  ’YDir’,  ’normal') 

cb_p=colorbar ; 

title (cb_p  ,  ’ p_ { i } ’ ) 

hold  on 

plot  (X,Y,  ’k’)  ; 
plot  (Y,X,  ’k’)  ; 
t  i  1 1  e  (  ’  Uniform  ’ ) 
xlabel ( ' Long .  ’ ) 
y label  (  '  Lat  .  ’ ) 


subplot (2 ,2 ,4) 

imagesc  ( [  -  4.5  ,4.5]  ,[  -  4.5  ,4.5]  ,  q_uni  ,  [0  ,  max(  q_uni  (:))]) 

set(gca,  ’YDir’,  ’normal’) 

cb_q=colorbar ; 

title  (cb_q  ,  ’  q_ { i  }  ’  ) 

hold  on 

plot  (X,Y,  ’k’)  ; 
plot  (Y,X,  ’k’)  ; 
t  i  1 1  e  (  ’  Uniform  ’ ) 
xlabel ( ’ Long .  ’  ) 
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y label  (  ’  Lat  .  ’ ) 

print  (  ’  Discrete_uni_uni  ’  ,  ’-dpng  ’  ) 

C.1.5  Discrete  Search  Function 

function  [  pd  ,  n  ,  imax  ,  Pd_cum2  ,k]  =  dis ere te_s e arc h2  (p  ,  q  ,  nT  , 
stop  ) 

°/cConducts  a  Discrete  Search  base  on  inputs  . 

%p=  Location  probability  density 
%q=Conditional  miss  probabilit  density 
°/cnT=max  number  of  glimpses 

%stop  =  Probability  of  Detection  to  stop  algorithm 
^Outputs 

%pd=updated  location  probability  density 
%n=search  matrix 
%imax  =  Last  box  searched 
%Pd_cum=  t  o  t  al  probabiliy  of  detection 
%k=number  of  glimpses 
if  -exist  (  '  stop  ’ ) 
stop=sum(p  ( : )  )  ; 
elseif  stop  >sum(p  ( : )  ) 
stop=sum(p  ( : )  )  ; 

end 

p0=l-sum(p  ( : )  )  ; 

n= zeros ( s ize (p ) ) ;%i  n  i  t i a  1  i  z e  look  counter 
Pd_cum  =  0; 

Pd_cum2=Pd_cum ; 

PC=P ; 
pc_0=p0 ; 

k  =  1 ; 

while  k<=nT  &&  Pd_cum2  ( end  )<  stop 
pd=pc  .*(1  -  q)  . *q . A n  ; 

[pmax  ,  imax  (k)  ]  =  max(pd  ( : )  )  ; 
pmax ; 
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[irow  ,  icol]  =  ind2sub(size(pd)  ,imax(k) ) 

30  n(irow  ,  icol  )=n(irow  ,  icol ) +  1; 

31  Pd_cum2(k+l)=Pd_cum2(k)+pmax ; 

32  k=k+l; 

33  end 

34  pd=pc  .*(1 -q)  .*q  .  An  ; 

35  end 


C.2  Bayesian 

C.2.1  Main  Code 

1  clc 

2  clear  all 

3  bayes_search_prob 

4  bayes_loc_probs2 

5  close  all 

6  %p-uniform 

7  bayes_pUni_qUni 
s  p0_uni_uni  =  pc_0 

s  bayes_pUni_qlongstrat 
io  p0_uni_ls=pe_0 
a  bayes_pUni_qbothstrat 

12  p0_uni_bs=pe_0 

13  %p-long  stratified 

14  bayes_plongstrat_qUni 

15  p0_ls_uni  =  pc_0 

16  bayes_plongstrat_qlongstrat 
n  p0_ls_ls  =pc_0 

is  bayes_plongstrat_qbothstrat 

19  p0_ls_bs  =pc_0 

20  %p-both  stratified 

21  bayes_pbothstrat_qUni 
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pO_bs_uni=pc_0 

23  bayes_pbothstrat_qlongstrat 

24  p0_bs_ls=pc_0 

25  bayes_pbothstrat_qbothstrat 

26  p0_bs_bs=pc_0 

27  %p-normal  center 

28  bayes_pnormwide_qUni 

29  p0_nw_uni  =  pc_0 

30  bay es_normwide_qlong strat 

31  p0_nw_ls=pc_0 

32  bayes_normwide_qbothstrat 

33  p0_nw_bs=pc_0 

34  %p- normal  side 

35  bayes_pnormside_qUni 

36  p0_ns_uni=pc_0 

37  bayes_normside_qlongstrat 

38  pO_ns_ls  =pc_0 

39  bayes_normside_qbothstrat 

40  p0_ns_bs  =  pc_0 

41  %p-normal  corner 

42  bayes_pnormc_qUni 

43  p0_nc_uni=pc_0 

44  b ay es_normc_qlong strat 

45  p0_nc_uni  =  pc_0 

46  b ay es_normc_qboth strat 

47  p0_nc_bs=pc_0 

48 

49  pd  =  [  Pd_uni_uni  ( end  )  ,  Pd_uni_ls  (  end  )  ,  Pd_uni_bs  (  end  )  ,Pd_ls_uni( 

end  )  , . .  . 

so  Pd_ls_ls  ( end  )  ,  Pd_ls_b s  (  end  )  ,  Pd_bs_uni  ( end  )  ,  Pd_b s_ls  (  end  ) 

, Pd_bs_bs ( end )  . .  . 

51  ,  Pd_nw_uni  ( end  )  ,  Pd_nw_ls  (  end )  ,  Pd_nw_bs  (  end )  ,  Pd_ns_uni  ( 

end  )  , .  .  . 
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52 


Pd_ns_ls ( end )  , Pd_ns_bs ( end  )  ,  Pd_nc_uni ( end )  ,  Pd_nc_ls ( end ) 
?  •  •  • 

53  Pd_nc_bs  ( end  )  ] 

54  m 

55 

56  close  all 

57  set(0,’defaulttextinterpreter’  , ’latex') 

58  [  ~  ,Pd_uni_uni  ,k,~  ,~]  =  bayes_search  (  P_uni  ,  q_uni  ,500  , 

•  95  ); 

59  [  ~  ,  Pd_uni_ls  ,k,~  ,~]  =  bayes_search  (  P_uni  ,  q_long_strat 

,500,  .95  ); 

60  [  ~  ,~  ,~  ,Pd_uni_bs  ,k,~  ,~]  =  bayes_search  (  P_uni  ,  q_both_strat 

,500,  .95  ); 

61  figure 

62  stepl  =10:  round(. 05*  length  (Pd_uni_uni))  :  length(Pd_uni_uni)  ; 

63  step2  =  10:  round  (.05*  length  (Pd_uni_ls))  :  length(Pd_uni_ls)  ; 

64  step3  =  10:round(.05*length(Pd_uni_bs))  :  length(Pd_uni_bs)  ; 

65  plot(stepl  ,Pd_uni_uni(stepl  )  ,  ’xk’  , step 2  ,Pd_uni_ls(step2)  ,  ’xr 

,  step  3  ,  Pd_uni_bs  (step3),’xg’) 

66  g  =  fittype  (’  a-exp  (b*x)  dependent  Prob  independent  ’,  ’x 

67 

68  [  fi  t_bb  1  ,  gof_bb  1  ]  =  fit  (stepl  ’  ,Pd_uni_uni(stepl  )  ’  ,  g  , 

StartPoint  ’  ,[0  ,0]  ,  ’  Upper  ’,[1,0]); 

69  [  fit_bb2  ,  gof_bb2  ]=fit(step2  ’  ,Pd_uni_ls(step2)  ’  ,g  ,  ’StartPoint 

’  ,[0,0]  ,  ’Upper’  ,[1  ,0])  ; 

70  [  fit_bb3  ,  gof_bb3  ]  =  fit  (step3  ’  ,  Pd_uni_bs  (step3)  ’  ,g,  ’  StartPoint 

’  ,[0,0]  ,  ’Upper’  ,[1  ,0])  ; 

71  hold  on 

72  plot  (  fit_bb  1  ,  ’ k  ’ ) 

73  plot  (  fit_bb2  ,  ’  r  ’ ) 

74  plot  (  f i t_b b 3  ,  ’  g  ’ ) 

75  xlabel(’n’) 
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y label ( ’$P_{D}$  ’ ) 

ylim  ( [0 , 1  ] ) 

title  (  ’  $P_ { D }  $  vs  Number  of  Glimpses  and  Exponential  Fitted 
Curves  ’ ) 

legend  (  ’  (p_{  uniform  }  ,q_{  uniform  })’,’( p_ {  uniform  }  ,q_{  strat  a 
})  ’  ,  ’  (p_{  uniform  }  ,q_{  strat  b  })  ’  .  .  . 

,  sprintf  (  '  %.5f-eA[t7o.5f*x]  ’  ,  f  i  t  _  b  b  1  .a,  f  i  t  _  b  b  1  .b)  ,  sprintf  ( 
’%.5f-eA{°/c.5f*x]  ’  ,  fit_bb2  .  a,fit_bb2  .b) 
sprintf  (  ,c7o.5f-eA{%.5f*x}  ’  ,fit_bb3  .a,  fit_bb3  .b)  ,  ’  location 
’  , ’ best ’ ) 

print  (  7  Bayes_Uni_Fit  ’  ,  ’-dpng  ’ ) 

m 

[  ~  ,  Pd_ls_uni  ,k,~  ,~]  =  bayes_search  (  P_long_strat  ,  q_uni 

,500,  .95  ); 

[  ~  ,~  ,~  ,  Pd_ls_ls  ,k,~  ,~]  =  bayes_search  (  P_long_strat  , 
q_long_strat  ,500 ,  .95  ); 

[  ~  ,  Pd_ls_bs  ,k,~  ,~]  =  bayes_search  (  P_long_strat  , 

q_both_strat  ,500 ,  .95  ); 

figure 

stepl=10:round(.05*length(Pd_ls_uni))  :  length(Pd_ls_uni)  ; 

step2  =  10:  round  (.05*  length  (Pd_ls_ls))  :  length(Pd_ls_ls)  ; 

step3  =  10:round(.05*length(Pd_ls_bs))  :  length(Pd_ls_bs)  ; 

plot(stepl  ,Pd_ls_uni(stepl)  ,’xk’  , step 2  ,Pd_ls_ls(step2)  ,’xr’  , 
step3  ,  Pd_ls_bs  (  step3  )  ,  ’  xg  ’ ) 

[  fi  t_bb4  ,  gof_bb4 ]=  fit  (stepl  ’  ,Pd_ls_uni(stepl  )  ’  ,g,  ’  StartPoint 
’  ,[1,0]  ,  ’Upper’  ,[1  ,0])  ; 

[  f  i  t_b  b  5  ,  gof_bb5  ]—  fit  (step2  ’  ,Pd_ls_ls(step2)  ’,g,  ’StartPoint  ’ 
,[1,0]  ,  ’Upper’  ,[1  ,0])  ; 

[  fit_bb6  ,  gof_bb6  ]=fit(step3  ’  ,Pd_ls_bs(step3)  ’,g,  ’StartPoint  ’ 
,[1  ,0]  ,  ’Upper’  ,[1  ,0])  ; 

hold  on 
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98  plot  (  fit_bb4  ,  ’k  ’ ) 

99  plot  (  f i t_b b 5  ,  ’  r  ’  ) 

100  plot  (  fit_bb6  ,  ’  g  ’ ) 

101  xlabel(’n’) 

102  ylabel  (  ’  $P_{D}$  ’ ) 

103  t i 1 1  e ( ’ $P_ {D} $  vs  Number  of  Glimpses  and  Exponential  Fitted 

Curves  ’ ) 

104  legend  (’ (p_{  strat  a  }  ,q_{  uniform  (p_{  s  tr  at  a},q_{strat  a 

})’,’(p_{strat  a},q_{strat  b})’... 

105  ,  sprintf(’%.5f-eA{%.5f*x}’  ,fit_bb4.a,fit_bb4.b),sprintf( 

’%.5f-eA{°/c.5f*x}  ’  ,  fit_bb5  .  a,fit_bb5  .b) 
ioo  sprintf(,%.5f-eA{c7c.5f*x}’  ,fit_bb6.  a,  fit_bb6.b),  ’location 

’  , ’ best ’ ) 

107  print  (  'Bayes_Long_Fit  ’  ,  ’-dpng  ’  ) 

108 

109  [  ~  , ~  ,  Pd_bs_uni  ,  k  ,  ~  ,  ~  ]  =  bayes_search  (  P_both_strat  ,  q_uni 

,500,  .95  ); 

no  [  ~  ,~  ,~  ,  Pd_bs_ls  ,k,~  ,~]  =  bayes_search  (  P_both_strat  , 
q_long_strat  ,500 ,  .95  ); 

in  [  ~  ,~  ,~  ,Pd_bs_bs  ,k,~  ,~]  =  bayes_search  (  P_both_strat  , 
q_both_strat  ,500 ,  .95  ); 
n2  figure 

n3  stepl=10:round(.05*length(Pd_bs_uni))  :  length(Pd_bs_uni)  ; 
n4  step2  =  10:round(.05*length(Pd_bs_ls))  :  length(Pd_bs_ls)  ; 
ns  step3  =  10:round(.05*length(  Pd_bs_bs  )):length(  Pd_bs_bs  )  ; 
lie  plot(stepl,  Pd_bs_uni  (stepl)  ,’xk’  , step  2  ,  Pd_bs_ls  (step2)  ,  ’  xr  ’  , 
step  3  ,  Pd_bs_bs  (step3)  ,’xg’) 

in  [  fit_bb7  ,  gof_bb7  ]=  fit  (stepl  ’,  Pd_bs_uni  (  stepl  )  ’  ,g,  ’  StartPoint 
’  ,[1  ,0]  ,  ’Upper’  ,[1  ,0])  ; 

ns  [  f  1  t_b  b  8  ,  gof_bb8  ]  =  fit  (  step2  ’  ,Pd_bs_ls(step2)  ’,g,  ’StartPoint  ’ 
,[1  ,0]  ,  ’Upper’  ,[1  ,0])  ; 

no  [  f  i  t_bb  9  ,  gof_bb9  ]=fit(step3’,  Pd_bs_bs  (step3)  ’  ,g,  ’  StartPoint  ’ 
,[1,0]  ,  ’Upper’  ,[1  ,0])  ; 
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hold  on 
m  plot  (  fit_bb7  ,  ’ k  ’  ) 

122  plot  (  fit_bb8  ,  ’  r  ’ ) 

123  plot  (  fit_bb9  ,  ’  g  ’ ) 

124  x label  (  ’n  ’  ) 

125  ylabel  (  ’  $P_{D}$  ’ ) 

126  title  (  ’  $P_ { D }  $  vs  Number  of  Glimpses  and  Exponential  Fitted 

Curves  ’ ) 

127  legend (’ (p_{  s tr at  b  },  q_  {  uniform  })’,’( p_  {  s  t  r  a  t  b},q_{strat  a 

})’,’( p_{strat  b  }  ,  q_  {  s  t  r  a  t  b  })’... 

128  ,  sprintf(’%.5f-eA[%.5f*x}’,fit_bb7.a,fit_bb7.b),sprintf( 

’°/c.5f-eA{°/c.5f*x}  ’  ,  fit_bb  8  .  a,fit_bb8  .b) 

129  sprintf(’%.5f-eA[%.5f*x}’,fit_bb9.  a,  fit_bb9.b),  ’location 

’  , ’ best ’  ) 

no  print  (  ’Bayes_Both_Fit  ’  ,  ’-dpng  ’ ) 

131 

132  [  ~  , ~  ,  Pd_nw_uni  ,  k  ,  ~  ,  ~  ]  =  bayes_search  ( P_norm_wide  ,  q_uni 

,500,  .95  ); 

133  [  ~  ,~  ,~  ,Pd_nw_ls  ,k,~  ,~]  =  bayes_search  ( P_norm_wide  , 

q_long_strat  ,500 ,  .95  ); 

134  [  ~  ,~  ,~  ,Pd_nw_bs  ,k,~  ,~]  =  bayes_search  ( P_norm_wide  , 

q_both_strat  ,500 ,  .95  ); 

135  figure 

136  stepl  =10:  round(. 05*  length  (  Pd_nw_uni  )):length(  Pd_nw_uni )  ; 

137  step2  =  10:round(.05*length(  Pd_nw_ls  )):length(  Pd_nw_ls  )  ; 
ns  step3  =  10:round(.05*length(  Pd_nw_bs  )):length(  Pd_nw_bs  )  ; 

139  plot(stepl,  Pd_nw_uni  (stepl)  ,’xk’  , step  2  ,  Pd_nw_ls  (step2)  ,  ’  xr  ’  , 
step  3  ,  Pd_nw_bs  (step3)  ,’xg’) 

mo  [  fit_bb  1  0  ,  gof_bb  10  ]=  f  i  t  (  step  1  ’  ,  Pd_nw_uni  (stepl)  ’  ,  g  , 
StartPoint  ’  ,[1  ,0]  ,  ’  Upper  ’,[1,0]); 
mi  [  fit_bb  1 1  ,  gof_bbl  1  ]=  f  it  (  step  2  ’  ,  Pd_nw_ls  (step2)  ’,g, 

StartPoint  ’  ,[1  ,0]  ,  ’  Upper  ’,[1,0]); 

W2  [  fit_bb  1  2  ,  gof_bb  12]  =  fit  (step  3  ’  ,  Pd_nw_bs  (step3)  ’,g, 
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StartPoint  ’  ,[1  ,0]  ,  ’  Upper  ’,[1,0]); 

143  hold  on 

144  plot  (  fit_bb  1  0  ,  ’k’) 

145  plot  (  fit_bb  11  ,  ’  r  ’ ) 
we  plot  (  fit_bb  1  2  ,  ’g’) 

M7  xlabel(’n’) 

148  ylabel  (  '  $P_{D}$  ’  ) 

149  title  (  ’  $P_ { D }  $  vs  Number  of  Glimpses  and  Exponential  Fitted 

Curves  ’ ) 

iso  legend  (’( p_{  norm  },  q_{  uniform  })’,’( p_{  norm  },  q_{  s  tr  at  a})’,’( 
p_  {  norm  }  ,  q_  {  s  t  r  a  t  b  } )  ’  .  .  . 

151  ,  sprintf(’%.5f-eA{c7o.5f*x}’  ,fit_bbl0.a,fit_bbl0.b)  , 

sprintf(’%.5f-eA[%.5f*x}’  ,  fit_bb  1  1  .a,fit_bbll  .  b ) 

152  sprintf(’^7o.5f-eA{c7o.5f*x}’,fit_bbl2.a,fit_bbl2.b),, 

location  ’  ,  ’  best  ’ ) 

153  print  (  ’  Bayes_Norm_Fit  ’  ,  ’-dpng  ’  ) 

154 

155  [  ~  ,~  ,~  ,  Pd_ns_uni  ,  k  , ~  , ~ ]  =  bayes_search  (  P_norm_side  ,  q_uni 

,500,  .95  ); 

156  [  ~  ,~  ,~  ,  Pd_ns_ls  ,k,~  ,~]  =  bayes_search  (  P_norm_side  , 

q_long_strat  ,500 ,  .95  ); 

157  [  ~  ,~  ,~  ,Pd_ns_bs  ,k,~  ,~]  =  bayes_search  (  P_norm_side  , 

q_both_strat  ,500 ,  .95  ); 

158  figure 

159  stepl  =10:  round(. 05*  length  (Pd_ns_uni))  :  length(Pd_ns_uni)  ; 
wo  step2  =  10:round(.05*length(Pd_ns_ls))  :  length(Pd_ns_ls)  ; 

161  step3  =  10:round(.05*length(  Pd_ns_bs  )):length(  Pd_ns_bs  )  ; 

162  plot  (  step  1  ,  Pd_ns_uni  (stepl)  ,’xk’  , step  2  ,Pd_ns_ls(step2)  ,  ’  xr  ’  , 

step  3  ,  Pd_ns_bs  (step3)  ,’xg’) 

163  [  fit_bb  1  3  ,  gof_bb  13]  =  fit  (step  1  ’,  Pd_ns_uni  (stepl)  ’  ,g  , 

StartPoint  ’  ,[1  ,0]  ,  ’  Upper  ’,[1,0]); 

164  [  fit_bb  1  4  ,  gof_bb  14  ]=  fit  (  step  2  ’  ,  Pd_ns_ls  (step2)  ’  ,g  , 

StartPoint  ’  ,[1  ,0]  ,  ’  Upper  ’,[1,0]); 
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165  [  fit_bb  1  5  ,  gof_bbl5  ]  =  fit  (step  3  ’  ,  Pd_ns_bs  (  step3  )  ’  ,g , 
StartPoint  ’  ,[1  ,0]  ,  ’  Upper  ’,[1,0]); 
we  hold  on 

167  plot  (  fit_bb  1  3  ,  ’k’) 

168  plot(fit_bbl4,’r’) 

169  plot  (  fit_bb  15  ,  ’  g  ’ ) 
no  x label  (  ’n  ’  ) 

Hi  ylabel  (  ’  $P_{D}$  ’ ) 

172  title  (  ’  $P_ { D }  $  vs  Number  of  Glimpses  and  Exponential  Fitted 
Curves  ’ ) 

m  legend  (’ (p_{norm  b  },  q_  {  uniform  })’,’( p_  {  norm  b},q_{strat  a}) 
, ’(p_{norm  b},q_{strat  b})’... 

174  ,  sprintf(’%.5f-eA{%.5f*x}’  ,fit_bbl3  .a,fit_bbl3  .b)  , 

sprintf  (  ’  %.5 f-e  A { %.5  f  *x }  ’  ,fit_bbl4  .a,  fit_bb  14  .b)  ,... 
ns  sprintf(’c7c.5f-eA{<7o.5f*x}’  ,  fit_bb  1  5  .a,fit_bbl5  .b)  ,’ 

location  '  ,  ’  best  ’  ) 

176  print  (  ’  B ayes_NormSide_Fit  ’  ,  ’-dpng  ’ ) 

177 

ns  [  ~  ,~  , ~  ,  Pd_nc_uni  ,  k  ,  ~  ,  ~  ]  =  bayes_search  (  P_norm_corner  ,  q_uni 
,500,  .95  ); 

no  [  ~  ,~  ,~  ,  Pd_nc_ls  ,  k  , ~  , ~ ]  =  bayes_search  ( P_norm_corner  , 
q_long_strat  ,500 ,  .95  ); 

iso  [  ~  ,~  ,~  ,Pd_nc_bs  ,k,~  ,~]  =  bayes_search  (  P_norm_corner  , 
q_both_strat  ,500 ,  .95  ); 

181  figure 

182  stepl=10:round(.05*length(Pd_nc_uni))  :  length(Pd_nc_uni)  ; 

183  step2  =  10:round(.05*length(Pd_nc_ls))  :  length(Pd_nc_ls)  ; 

184  step3  =  10:round(.05*length(  Pd_nc_bs  )):length(  Pd_nc_bs  )  ; 

185  plot(stepl,  Pd_nc_uni  (stepl)  ,’xk’  , step  2  ,Pd_nc_ls(step2)  ,’xr’ 

step  3  ,  Pd_nc_bs  (step3)  ,’xg’) 

186  [  fit_bb  1  6  ,  gof_bb  16]  =  fit  (step  1  ’,  Pd_nc_uni  (stepl)  ’,g, 

StartPoint  ’  ,[1  ,0]  ,  ’  Upper  ’,[1,0]); 
is?  [  fit_bb  1  7  ,  gof_bb  17  ]  =  fit  (step  2  ’  ,Pd_nc_ls(step2)  ’,g, 
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StartPoint  ’  ,[1  ,0]  ,  ’  Upper  ’,[1,0]); 

[  fit_bb  1  8  ,  gof_bb  18]  =  fit  (step  3  ’,  Pd_nc_bs  (step3)  ’  ,g  , 

StartPoint  ’  ,[1  ,0]  ,  ’  Upper  ’,[1,0]); 
hold  on 

plot  (  fit_bb  1 6  ,  ’k’ ) 
plot  (  fit_bb  17  ,  ’  r  ’ ) 
plot  ( fit_bb  1  8  ,  ’g  ’ ) 
xlabel(’n’) 
y  label ( ’$P_[D}$ ’ ) 

t  i  1 1  e  (  ’  $P_ [D }  $  vs  Number  of  Glimpses  and  Exponential  Fitted 
Curves  ’ ) 

legend  (’ (p_{  norm  c  },  q_{  uniform  })’,’( p_{  norm  c],q_[strat  a})’ 
,’(p_{norm  c],q_[strat  b})’... 

,  sprintf  (  ’  %.5  f-e  A{  %.5  f  *x }  ’  ,fit_bbl6  .a,  fit_bb  1  6  .b)  , 

sprintf(  ’  %.5 f-e  A{ 90.5  f  *x }  ’  ,  fit_bb  1  7  .a,  fit_bb  1  7  .b)  ,... 
sprintf  (  ’  %.5 f-e  A { 9o.5  f  *x }  ’  ,  fit_bb  1  8  .a,  fit_bb  1  8  .b)  ,  ' 
location  ’  ,  ’  best  ’  ) 

print  (  ’  B ayes_NormCorner_Fit  ’  ,  ’-dpng  ’ ) 
format  shortG 

b  =  [fit_bbl  .b,  fit_bb2  .b,  fit_bb3  .b,  fit_bb4  .b,  f i t _ b b 5  .b,  fit_bb6 
.  b  ,  fit_bb7  .  b  ,  .  .  . 

fit_bb 8  . b  ,  fit_bb9  .b,  fit_bb  1 0  . b  ,  f it_bb  1 1  .b,  fit_bb  1  2  .b, 
fit_bb  1  3  .  b  ,  fit_bb  1  4  .b  ,  .  .  . 
fit_bb  1 5  . b  ,  fit_bb  1  6  .b,  fit_bb  1 7  .b,  fit_bb  1  8  .b] 
r2  =[ gof_bb 1  .rsquare  , gof_bb2 . r square  , gof_bb3 . r square  , gof_bb4 . 
rsquare  ,  .  .  . 

gof_bb5  .  rsquare  ,  gof_bb6  .  rsquare  ,  gof_bb7  .  rsquare  ,  gof_bb8  . 
rsquare  ,  .  .  . 

gof_bb9  .  rsquare  ,  gof_bbl0  .  rsquare  ,  gof_bbl  1  .  rsquare  , 
gof_bbl2  .  rsquare  ,  .  .  . 

gof_bbl3  .  rsquare  ,  gof_bbl4  .  rsquare  ,  gof_bbl5  .  rsquare  , 
gof_bbl6  .  rsquare  ,  .  .  . 
gof_bbl7  .  rsquare  ,  gof_bbl8  .  rsquare  ] 
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[  ~  ,Pd_uni_bs  ,k,~  ,~]  =  bayes_search  (  P_uni  ,  q_both_strat 
,500,  1); 

[  ~  ,  Pd_uni_bs_r  ,k]  =  bay es_search_random  (  P_uni  , 

q_both_strat  ,500  ,  1  ); 
figure 

plot  (1:  length(Pd_uni_bs)  ,  Pd_uni_bs  ,’xr’  ,  1 :  length  (Pd_uni_bs_r 
)  ,  Pd_uni_bs_r ,  ’ r ’ ) 

Wc 

for  i  =  1 : 9 

p0  =  .  1  *  i  ; 

P_uni  =  ones  (10,10)  ; 

P_uni  =(l-p0)*P_uni/ sum(  P_uni  ( : )  )  ; 

[  ~  ,~  ,~  ,  Pd_uni_bs_p0  ( i  , : )  ,k,~  ,~]  =  bayes_search  (  P_uni  , 
q_both_strat  ,500  ,  1); 

end 

figure 
for  i  =  1 : 9 
hold  on 

plot  (1:25:500  ,  Pd_uni_bs_p0  (i  ,1:25:500)) 

end 


C.2.2  Sample  code  for  the  Uniform/Uniform  case 

[  pc  ,  n  ,  imax  ,  Pd_uni_uni  ,k  ,  pc_0  ,  k_check  ]  =  bayes_search  (  P_uni 
,  q_uni  ,50  ,  .95  )  ; 
figure 

subplot  (2,2,1) 

imagesc  ( [  -  4. 5, 4. 5], [-4. 5, 4. 5],  pc, [0  ,  max  (pc  (:))]) 

set(gca,  ’YDir’,  ’normal’) 

cb_p=colorbar ; 

title (cb_p  ,  ’p_{ i } ’ ) 

hold  on 
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plot  (X,Y,  ’k’)  ; 

plot  (Y,X,  ’k’)  ; 

title  (  ’  Remainder  Prob  .  ’ ) 

xlabel ( " Long .  ’ ) 

y label  (  ’  Lat  .  ’  ) 


subplot(2,2,2) 

imagesc  ( [  -  4.5  ,4.5]  ,[  -  4.5  ,4.5]  ,n,[0  ,  max(n  (:))]) 

set(gca,  ’YDir’,  ’normal’) 

cb_p=colorbar ; 

title  ( cb_p  ,  ’N’  ) 

hold  on 

plot  (X,Y,  ’k’)  ; 
plot  (Y,X,  ’k’)  ; 
t  i  1 1  e  (’ Actual  Search’) 
xlabel ( ’ Long .  ’  ) 
y label  (  ’  Lat  .  ’  ) 
subplot (2 ,2 ,3) 

imagesc  ( [  -  4.5  ,4.5]  ,[  -  4.5  ,4.5]  ,  P_uni  ,  [0  ,  max(  P_uni  ( : )  )  ] ) 

set(gca,  ’YDir’,  ’normal’) 

cb_p=colorbar ; 

title (cb_p  ,  ’ p_ { i } ’ ) 

hold  on 

plot  (X,Y,  ’k’)  ; 
plot  (Y,X,  ’k’)  ; 
t  i  1 1  e  (  ’  Uniform  ’ ) 
xlabel ( ' Long .  ’ ) 
y label  (  ’  Lat  .  ’  ) 


subplot  (2,2  ,4) 

imagesc([  -  4.5  ,4.5]  ,[  -  4.5  ,4.5]  ,  q_uni  ,  [0  ,  max(  q_uni  (:))]) 
set(gca,  ’YDir’,  ’normal’) 
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cb_q=colorbar ; 
title  (cb_q  ,  ’q_{ i } ’ ) 
hold  on 

plot  (X,Y,  ’k’)  ; 
plot  (Y,X,  ’k’); 
title  (  ’  Uniform  ’ ) 
xlabel ( ’ Long .  ’  ) 
y label  (  ’  Lat  .  ’  ) 

print  (  ’Bayes_uni_uni  ’  ,  ’-dpng  ’ ) 

C.2.3  Bayes  Search  Function 

function  [  pd  ,  n  ,  imax  ,  Pd_cum  ,k  ,  pc_0  ,  k_check  ]  =  bayes_search  ( 
p,q,nT,  stop  ) 

^Conducts  a  Bayesian  Search  base  on  inputs  . 

%p=  Location  probability  density 
%q=Conditional  miss  probabilit  density 
%nT=max  number  of  glimpses 

%stop  =  Probability  of  Detection  to  stop  algorithm 
%Outputs 

%pd=  updated  location  probability  density 

%n=search  matrix 

%imax  =  Last  box  searched 

%k=number  of  glimpses 

%pc_0=  f  in  al  probability  tgt  is  not  in  search  area 
%k_check=game  theory  detection  glimpse 
if  -exist  (  '  stop  ’ ) 
stop=sum(p  ( : )  )  ; 
elseif  stop  >sum(p  ( : )  ) 
stop=sum(p  ( : )  )  ; 

end 

check  =0; 

n= zeros ( s ize (p ) ) ;%i  n  i  t i a  1  i  z e  look  counter 
Pd_cum  =  0; 
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pc_0=l-sum(p  ( : )  )  ; 

PC=P ; 

k=  l ; 

tgt_pos  =r  and  sample  (1:101  ,1  ,true  ,[p(:)  ;  pc_0  ] )  ; 

[  tgt_row  ,  tgt_col]=ind2sub(  size(p)  ,tgt_pos  )  ; 
while  k<=nT  &&  Pd_cum  ( end  )<  stop 
pd=pc  .*(  1  -q)  ; 

[pmax  ,  imax  (k)  ]  =  max(pd  ( : )  )  ; 

[irow  ,  icol]  =  ind2sub(size(pd)  ,imax(k) )  ; 
if  check==0 

if  irow==tgt_row  &&  ic o  1  ==  t g t_co  1 

check  =  binornd(l,l  -  q(  irow  ,  icol  ) )  ; 
if  check==l 

k_check=k ; 

end 

end 

end 

pc  (irow  ,  icol  )=pc(irow  ,  icol  )  .  *q(irow  ,  icol )  ; 

pc=pc/(l  -pmax)  ; 

pc_0=pc_0/(l  -pmax)  ; 

n(irow  ,  icol  )=n(irow  ,  icol )  +  1; 

Pd_cum(k  + 1  )=Pd_cum  (k)  +(1  -Pd_cum(k) )  *pmax  ; 
k=k+  1 ; 

end 

if  check==0 

k_check=NaN ; 

end 

pd=pc ; 
end 


C.3  Continuous 

Code  is  as  in  [6]  with  the  exception  that  R\  =  Roe~vl^v  and  =  Roe~V2^v 
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